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mten Aufsätze der EncyklopKdie ebenso wie die Eeferate 
^.jberichten der Deutschen Mathematiker -Vereinigung be- 

in dieaem Sinne in knapper, füi- eine rasche Orientierung 
bestimmter Form den gegenwartigen Inhalt einer Diaciplin an gesicherten 
Besnltaten zn geben, wie auch durch sorgfältige Litteraturangaben die 
historische Entwickelung der Methoden darzulegen. Darüber hinaus aber 
mufs auf eine eingehende, mit Beweisen versehene Darstellung, wie sie . 
zum selbständigen, von umfangreichea Quellenstudien unabhängigen Ein- 
dringen in die Disciplin erforderlich ist, auch bei den breiter angelegten 
Beferateu der Deutscheu Mathematiker -Vereinigung, in welcher haupt- 
sächlich das historische und teilweise auch das kritische Element zur 
Geltuug kommt, Terzichtet werden. Eine solche ausfahrliche Darlegung, 
die sich mehr in dem Charakter eines auf geschichtlichen uud littera- 
rischeu Studien gegründeten Lehrbuches bewegt und neben den rein 
wissenschaftlichen auch pädagogische Interessen berücksichtigt, erscheint 
aber bei der raschen Entwickelung und dem umfang des zu einem 
grofsen Teil nur in Monographien niedergelegten Stoffes durchaus wichtig, 
zumal, im Vergleiche z. B. mit Frankreich, bei »ms in Deutschland die 
mathematische tiitteratur an Lehrbüchern über spezielle Ciebiete der 
mathematiBchen Forschung nicht aUzu reich ist. 

Die Verlt^buchhandlung B. G. Teubner giebt sich der HofEuung 
hin, dafs sich recht zahlreiche Mathematiker, Physiker und Astronomen, 
Geodäten und Techniker, sowohl des In- als des Auslandes, in deren 
Forschungsgebieten der;irtige Arbeiten erwünscht sind, zur Mitarbeiter- 
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Schaft an dem UntemeluiLeii entscbliefsen möchten. Besonders nahe liegt 
die Beteiligung den Herren Mitarbeitern an der Encyklopädie der Mathe- 
matischen Wissenschaften. Die mnfangreichen litterarischen und speziell 
fachlichen Studien, welche f£Lr die Bearbeitung von Absphnitten der 
Encyklopädie vorzunehmen waren, konnten in dem notwendig eng be- 
grenzten Bahmen nicht vollständig niedergelegt werden. Hier aber, bei. 
den Werken der gegenwärtigen Sammlung, ist die Möglichkeit gegeben, 
den Stoff freier zu gestalten und die individuelle Auffassung und Richtung 
des einzelnen Bearbeiters in höherem Mafse zur Geltung zu bringen. 
Doch ist, wie gesagt, jede Arbeit, die sich dem Plane der Sammlung 
einfügen läTst, im gleichen MäXse willkommen. 

Bisher haben die folgenden Gelehrten ihre geschätzte Mitwirkung 
zugesagt, während erfreuHcherweise stetig neue Anerbieten zur Mitarbeit 
an der Sammlung einlaufen, worüber in meinen „Mitteilungen*^ fortlaiifend 
belichtet werden wird (die bereits erschienenen Bände sind mit zwei **, 
die unter der Presse befindlichen mit einem * bezeichnet): 
^^F« Baohmanp, niedere Zahlentheorie. (Band X der Sammlung.) 

M. Böoher, über die reellen Lösimgen der gewöhnlichen linearen 
Differentialgleichungen zweiter Ordnung. 

G. BohlTnann, Yersicl^erungsmathematik. 

G. H. Bryan, Lehrbuch der Thermodynamik. 

G-. Castelnuovo und F. Enriques, TheoHe der algebraischen Flächen. 
'^E. Czuber, Wahrscheinlichkeitsrechnung und ihre Anwendung auf Fehler^ 

ausgleichung, Statistik und Lebensversicherung. (Band IX.) 
^X. E, Diokson, Linear Groups with an exposition of the Galois Field 

theory. [in engUsoher Sprache.] (Band VI.) 

F. Dingeldey, Kegelschnitte und Kegelschnittsystemew 

F. Dingeldey, Sammlung von Aufgaben zur Anwendung der Differential- 

imd Litegralrechnung. 

G. Eneström (in Verbindung mit andern Gelehrten), Handbuch der 

Geschichte der Mathematik. 

F. Enriques, Prinzipien der Geometrie. 

Ph. Fortwängler, die Mechanik der einfachsten physikalischen Apparate 

imd Versuchsanordnungen. 
**A, Gleichen, Optische Abbildungslehre u. Theorie der optischen 

Instrumente. (Band VIII.) 
M. Grübler, Lehrbuch der hydraulischen Motoren. ^ 

J. Harkness, elliptische Funktionen. 
L. Henneberg, Lehrbuch der graphischen Statik. 
K. Heun, die kinetischen Probleme der modernen Maschinenlehre. 

G. Jung, Geometrie der Massen. 
G. Kohn, rationale Kurven. 

'^A. ExaEer, Handbuch der Lehre von den Thetafunktionen. (Band YII*) 
H. Lamb, Akustik. 
B. V. Iiilienthal, Differentialgeometrie. 



A. Loewy, Vorlesungen über die Theorie der linearen Substitntionsgruppen. 
'^G. Loria, spezielle, algebraisclie nn^ transcendente Kurven der Ebene. 
Theorie und Geschichte. (Band V.) 
A. E. H. Iiovö, Lehrbuch der Hydrodynamik. 

A. E. H. Lov6, Lehrbuch der Elasticität. 

B« Mehmke, über graphisches Rechnen und über Rechenmaschine 

sowie Über numerisches Eechnen. 
W. Meyerhofer, die mathematischen Grundlage^ der Chemie. 

1. Netto, Lehrbuch der Combinatorik. (Band VII.) 
W. F. Osgood, aligemeine Funktionentheorie. 

E. Ovazza, aus dem Grebiete der Mechanik. 

**E. Fasoal,. Determinanten. Theorie und Anwendungen. (Band HL) 

S. Finoherle, Funktional -Gleichungen imd -Operationen. 

Fr. il^ookels, Kzystalloptik. 

A* Pringsheim, Vorlesungen über Zahlen- und Funktionenlehre. ^ (Ele- 
mentare Theorie der unendlichen Algorithmen imd der analytischen 
Funktionen einer komplexen Veränderlichen.) Bd. L Zahlenlehre. 
Bd. n. Funktionenlehre. (Band I.) 

O. Segre, Vorlesungen über algebraische Geometrie, mit besonderer 
Berücksichtigung der mehrdimensionalen Bäume. 

D. Seliwanoff, Differenzenrechnung. ^ 
M. Simon, Elementargeometrie. 

F« St&okel, Lehrbuch der allgemeinen Dynamik. 

F. Stäokel, Differentialgeometrie höherer Mannigfaltigkeiten. 
O. Staude, Flächen und Flächensysteme zweiter Ordnung. 

**0. StoLs und J. A. Gmeiner, theoretische Arüihmetik. (Band IV.) 

B. Sturm, Theorie der geometrischen Verwandtschaften. 
B. Sturm, die kubische Baumkurve. 

H« E. Timerdiug, Theorie der Streckensysteme und Schrauben. 
K. Th. Vahlen, Geschichte des .Fimdamentalsatzes der Algebra. 

E. Th.. Vahlen, Geschichte des Sturmschen Satzes. 
A. Voss, Prinzipien der rationellen Mechanik. 

A, Voss, Abbildimg und Abwicklung der krummen Flächen. 
J. G. Wallentüi, Lehrbuch der theoretischen Elektrizität. ^ 
**B, V, Weber, Vorles^ngen über das Pfaffsche Problem und die- Theorie 
der partiellen Differentialgleichimgen 1. Ordnung. (Band 11.) 
*A. G. Webster, the Dynamics of Pariicles, of rigid, elastic, and fluid Bodi^s 
beingLecturesonMathematicalPhysics. [in engUBoher Sprache.] (Band XI.) 
A. Wiman, endliche Gruppen linearer Transformationen. 
w, Wirtinger, algebraische Funktionen und ihre Integrale. , 
W. Wirtinger, partielle Differentialgleichungen. 
H. G. ZexLthen, die abzählenden Methoden der Geometrie. 

]t^ Mitteilungen über weitere Bände werden baldigst fplgen. 
Leipzig, Poststrafse 3; *> rs. m^««i^««^«. ' 

Oktober 1902. B, G. TGUDner. 
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V. 



Vorwort. 



— 

-Die älteren Werke über Dififerenzenreclmung besitzen nicht die 
Strenge^ die nnsrer Zeit eigentümlich ist. In dieser Beziehung ist 
befriedigend das Bach von A. A. Markoff (Differenzenrechnung^ 
St. Petersburg 1889—1891; deutsch von Th. Friesendorf und 
E. Prümm. Leipzigs bei B. G. Teubner. 1896). Ich habe versucht 
einige Verbesserungen zu machen. Vielleicht ist es mir gelungen 
u^ einige Beweise zu vereinfachen und den Stoff systematischer an- 

zuordnen. Um die Arbeit den Studierenden zu erleichtern^ habe ich 
speziellere Probleme beiseite gelassen und mich hauptsächlich auf 
elementare Fragen beschränkt. 

Worin nach meiner Meinung das WesentUche in der Differenzen- 
rechnung besteht^ habe ich in einem Referat dargestellt. (Encyklo- 
^ .o pädie der mathematischen Wissenschaften, Bd. I, p.918 — 937.) Dort 

findet man hinreichende Literaturangaben. In kleinem Baume konnte 
ich nur Sätze aussprechen, die Beweise sollte der Leser selbst suchen. 
Um ihm die erforderliche Mühe und Zeit zu ersparen, folge ich 
dem Vorschlag des Herrn B. G. Teubner und übergebe der Öffentlich- 
keit dieses Lehrbuch. 

St. Petersburg, den 6. Dezember 1903. 



D. Seliwanoff. 
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Differenzen, 



Erstes Kapitel. 

Hauptsätze über Differenzen. 

1. Definitionen. Wenn eine Funktion f{x) gegeben ist, so 
heißt der Ausdruck f{x + K) —f{x) Differenz von f{x)f in Zeichen: 

fix + A) - fix) = t.f{x). 

Wenn man diesen Ausdruck in eine andre Form bringt, so 
entsteht die Gleichung 

fix + A) - fix) = q>{x) 

oder 

t^f{x)^q)ix). 

Hier kann x beliebige Werte annehmen, A soll aber konstant 
bleiben. 

Ist z. B. fix) = x^, A = 1, so hat man 

/^(0) = 0, m) = l, f(2) = 8, f(3) = 27, /•(4) = 64,... 

A/'CO) = 1, £^f(^ = 7, Af(^= 19, ^m= 37, . . . 

Die Differenz von A/"(a5) heißt zweite Differenz von /"(«) oder 
Differenz zweiter Ordnung: 

AY(a:) = A/(a; + A) - A/^(«). 
Entsprechend ist 
AY(a;) = AY(a; + A) - AY(a;), AY(«) = AYC« + Ä) - AY(a;), . . . 
In dem Beispiel ist 

AY (0) = 6, AY(1) = 12, AY(2) = 18, . . . 
AY(0) = 6, AY(1) = 6,... 

AY(0) = 0,... 

Seliwanoff, Differenzenreclmting. 1 
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Diese Besnltate lassen sich in folgender Tabelle zusammenstellen: 
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2. Dllforensen einflMher Fnnktioiien. unter einfachen 
Funktionen versteht man diejenigen, deren Differenz direkt bestimmt 
werden kann. Die Funktion heißt zusammengesetzt, wenn ihre 
Differenz dorch die Differenzen andrer Funktionen ausgedrückt 
werden kann. 

In der Differentialrechnung ist 

Die Differenz von re" hat keine einfache Form, weil 

Aa!" = nÄa?— > _^ n(n-l) ^g^_a _^ ^ ^„ 

1 * « 

ist. Die einfache Funktion^ die der Potenz entspricht, ist in der 
Differenzenrechnung die folgende 



f{x) = x{^ — }i){x — 2A) . . . (a; — n — lA). 
Die direkte Rechnung ergibt 



fix + A) — fix) = (a; + }i)xix — A) . . . (a: — n — 2A) 

— ir(a: — A)(ic — 2A). . . (a; — w — lA). 

Um diesen Ausdruck zu transformieren, setzen wir Tor die 
Klammem den gemeinsamen Faktor ->, 

a?(a;-A)(a;-2A)...(a?-;r^A)(x^/» X-^^^-^--/) 

und machen Reduktionen in den Klammem. 
Auf diese Weise findet man die Formel 



b\xix — h)(x — 2A) . . .(x — n — lA)] 



= nhx(x — A)(a; — 2A) . . . (a? — n — 2A). 

Spezielle Fälle dieser Formel sind 

Aa? = h, A[a;(a? — A)] = 2Aaj, 

A[a;(a; — A)(a; — 2A)] = Shx{x — A), . . • 



Erstes Kapitel. Hauptsätze über Differenzen. 
Etwas einfacher ist die Formel 



^x(x — h)(x—2h). ..(rc — n— lÄ) , x{x — h){x — 2h)...(X'- n — 2h) 

12. 3. ..n "~ 1.2.3.. .(n-1) ' 

die ganz analog bewiesen wird. 

Wir werden noch folgende einfache Funktionen angeben: 



x{x + h){x + 2h)...{x + n'-lh) 

Hier ist m eine Eonstante. 

Die direkte Rechnung ergibt 

1 



} m^j siniz?, cosic. 



a;(a? + Ä)(aj + 2Ä)...(a; + n-.lÄ) (a? + Ä)(a; + 2Ä)(a; + 3Ä) . . . (aj + nÄ) 

1 > 

a? (a; + Ä) (a? + 2 Ä) . . . (a? + n — 1 Ä) 

= \ f-^ 11 

(aj + Ä)(« + 2Ä) . . . (a! + n-lÄ)L« + »»A «J 

** ' aj(a; +Ä)(a? + 2Ä) . . . (aj + nÄ)* 
Spezielle Falle dieser Formel sind 



X xipi^lC) x{pii-\-lC) a5(aj4-Ä)(a5-f 2Ä) 

Etwas symmetrischer ist die Formel 

. 1.2.8...(n— 1) , 1-2 3.. .n 

a;(aj + Ä)(aj + 2Ä)...(a; + n^Ä)"" 'a:(aj + Ä)(a; + 2Ä) . . . (a; + nÄ)' 

Eine sehr einfache Funktion ist m% weil 

A w* = w*+* ~ m* = w*(w* — 1). 

Um die Differenz dieser Funktion zu bilden^ genügt es^ m^ mit 
dem konstanten Faktor m^ — 1 zu multiplizieren. 

Nach den bekannten Formeln der Trigonometrie ist 

sin (aj + Ä) — sin ic = 2 sin— cos (a: + ^^j = 2 sin — sin (a? + ^^ + ^)' 

cos(a; + A) — cosa;=» — 2sin-2 sii^ (^ +-2) ~^^^¥^^®(^ + 2^ + Y/ 
Also wird 

Asina?= 2 sin— cos (a; + 0^) =* 2 sin — sin \x -\ — J^)' 

Acosa; = — 2sin — sin (a? + ^) == 2sin — cosfa? H — y—)* 



1 



« 
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Die wiederholte Anwendung der angegebenen Formeln liefert 
uns folgende Ausdrücke für Differenzen höherer Ordnung 



^ 1.2.8...n * 1.2.8...(n-m) ^»»<«; 

. n X{x — h) (« — 2Ä). . .(x — n — lh) ,n 

^ 1.2.8...n —"■ ' 



m 



1.2-3. ..(n—1) , ,xn» 1-2.3. ..(w + w—1) 



x{x + Ä) (a? -f- 2Ä) . . . (a; + n — lÄ) x(x -f- ä) . . . (rc + n + m — lÄ) 

w = m (^w — 1; , 
A'*sina; = (2 sin^j • sin (o? + « • -y^)' 

A^cosic = (2 sin— j . cos (a? + w • — |— ) • 

Das sind die wichtigsten Formeln ^ die für das Folgende aus- 
reichen. 

3. Differenzen zusammengesetzter Funktionen. Wenn 

TTobei C eine Konstante ist, so ist 

f{x + Ä) - f{x) = [<jp(a; + Ä) + CQ - [y(a;) + C] = ,,(3; + Ä) - <p(a;), 

also Af{x)^Aq){x). 

Bei der Berechnung der Differenz kann also der konstante Sum- 
mand weggelassen werden. 

ist, so ist 

Af{x) = Cq>{x + h)- Cq){x) = Clq>(x + h)- (p{x)'] 

oder A [Cy («)] = A 9) (a;). 

Der konstante Faktor kann also Tor das Zeichen A gesetzt 
werden. 

Ist die gegebene Funktion f(x) Samme von mehreren Funktionen: 

fix) = f,(x) + f,{x) + . . . + Ux), 
so findet man 

^f(x)=[f,(x+h)+f,ix+h)+- ..+f^(x+h)]-\f,(x)+f,(x)+- ■ ■+Ux)] 
={fii^+1^)-fi{'»y] + [f»(x + h)-f,(xy] + -- •+Mx + h) -Uix)]. 
Es besteht also die Formel 
^[/i(«) + /»(«') + • • • + UixJ] = A/iC») + Af.ix) + ■-' + AUx). 
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Die bewiesenen Sätze lassen sich in einer Formel vereinigen. 
Es ist 

wobei C, Gl, G^, . , ,yGfi konstant sind. 

4. Differenzen ganzer Funktionen. Es sei 

Es ist bequemer, diese Funktion zuerst in die Form 

r./ \ A \ A oc — a , . (X'-a)(x—a—h) , » (x—a)(x—a-~h)(x—a—2h) . 

fix) = A, + A,-^ +^^ 'j^ '- +^,^ 1.2.3 + •■■ 



. j, {:f — a){x — a — h) . . . (x — a — n — Ih) 
••■+^ 1.2.3...n 

ZU bringen. Zur Bestimmung der Koeffizienten bemerken wir, daß 
bei der Division von f(x) durch x — a der Rest Äq ist und der 

Quotient 

/./x ^ , s x — a — h , . (x — a — h)(x — a — 2h) . 

fM=Ä, + Ä,-^^j-+A^'^ 1^ -+■■■ 



.j (x — a — h)...{x — a — n-~lh) 
""^^ 1.2.3...W 

Es ist noch zu bemerken, daß fi(x) die Form hat 

Dividiert man f^(x) durch x — a — h, so findet man den Rest Ä^ 
und den Quotient 



hW-AlT^+A 1.2.3 +--- + ^- 1.2.3...n 

Das erste Glied von f^ix) ist p^a:""-*. 

Schließlich ist 

r / \ A 1 i A a;— a— n — lÄ , 

/•«-i(a;)=A-i- i.g (^_i) +A 1.2. 3. ..n -P«^ + «- 

Hier ist g eine Zahl, die sich vollständig bestimmen läßt. 

Durch Division von fn—i(x) durch x — a — n — 1 findet man 

^ ^ A 

den Rest :r— — —, — tt und den Quotient ^ ^ ^ — Dieser Quotient ist 

l-2...(w — 1) l'2...n 

aber ^q. Es ist also 

1.2.3...n"~^o* 

Nachdem die Koeffizienten -^q, -4i, A^, . , ,, An bestimmt sind, 
erhält man mit Anwendung der Formeln (Nr. 2): 
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Afix) ^h[Ä, + A^-^+-' + A^ l.2.3...(n~l) -y 

Af{x) = h^[A^ + A^-j-+'- + A^ 1.2.3...(n-2) ^J' 

Ay{x)^'hr'Än 

oder 

Es bestellt also der Satz: 

Die erste Differenz einer ganzen Funktion w*®^ Grades 
ist vom (w — 1)*®^ Grade und die n^ Differenz ist konstant. 

Für a; = a ist 

/•(a)=A, Af{a)^h'A> A?f{a) = h^.A,,... 

Arf(a) = 'hr'An. 

Wir sehen^ daß die Koeffizienten A^, Ä^, , . ., An sich durch 
Differenzen ausdrücken lassen: 

Diese Werte der Differenzen bestimmen sich durch die Werte 
der Funktion f{x): 

f(a), fia + h), f(a + 2h),...J(a + nh). 

Da aber f(x) eine beliebige ganze Funktion n*^^ Grades ist, so 
besteht die Formel 

t(x)^f {a) + — j^ + ^ ji- + . • . 



"i 1.2-8.. . w ' ^i 

für jede Funktion von der Form 

wobei einige Koeffizienten gleich Null sein können. 
Das ist die Newtonsche InterpolationsformeL 
Es sei z. B. 

a = 0; Ä = 1, f{x) = a;*. 

Wir haben schon in Nr. 1 gesehen, daß 

/•(0) = 0, A/(0) = 1, AY(0) = 6, AY(0) = 6. 
Polglich ist 

x^ = x-\- Zx(x - 1) + x{x — !)(« — 2). 
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Als zweites Beispiel nehmen wir 

a-ü, Ä-1, t(x)^ i.2.3...n 

Hier ist y konstant vorausgesetzt. 

Nach der bewiesenen Formel (Nr. 2) ist 

'^^ 1.2.3...(n--w) 

Daraus folgte daß 

(x + y)(x + y-l)...{x + y-n+l) ^ y(y-l)...(y-n+l) 
1-2- 3 . ..n 1.2-3.. .n 

, X y(y-l),..(y-n + 2) a;(a;~l) y(y-l)...(y-n + 3) 



+ 



1 1.2.3...(n-l) ' 1-2 1-2. ..(n-2) 

35(0?— 1) .. .(o; — n+l) 



l-2-3...n 

Diese Relation bestellt für beliebige Werte von x und y; sie ist von 
Cauchy benutzt worden, um den binomischen Lehrsatz auf beliebige 
reelle Exponenten auszudehnen.^ 

6. Differens des Produktes sweier Funktionen. Wenn 
man die Differenzen zweier Funktionen q)(x) und ^(x) kennt, so ist 
es leicht, die Differenz ihres Produktes zu bilden. 

Wenn f(x) == q>{x)' ^{x) ist, so ist 

^f{x) =:q)(x + h)' tl;(x + h) — (p{x) • f{x) 

= IvQ^) + A9?(a;)]^(ir + h) - tp{x) • ^{x) 
= q>{x)\:^{x + A) - ^{xy] + i){x + A) . A9?(a?). 
Man findet also die Formel 

£^\sp{x) • ^(o?)] = (p{x) • A^(ir) + ^{x + h) 'Aq){x). 
Z. B. für A = 1 ist 

6. Differensen von sin u und cos u. Wir wollen voraussetzen, 
daß u eine Funktion von x ist, deren Differenz bekannt ist. Es ist 
Asinte und Acost^ zu berechnen. 

Wenn man x durch x + h ersetzt, so geht w in w + £>^u über. 
Deswegen ist 

smM=: sm(M + Aw) — smw = 2sin-r- • cosfw +-0")' 
Acosw = cos(w + Au) — cosw = — 2sin-2^ • sin (u +-y)* 



Cauchy, Conrs d'analyse, p. 159. Paris. 1821. 
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Damit besclili'eßen wir die Sätze über Differenzen zusammen- 

• 

gesetzter Funktionen. 

7. Berechnung Aquidlstanter Werte einer ganzen Funk- 
tion. Es sei f{x) eine ganze Funktion n^^ Grades, deren Werte 

ausgerechnet seien. 

Es sollen die Werte 



f(a + n+lh), f(a + n + 2h), /"(a + w + 3ä), . . . 
und 

f(a-h), f{a-2h), f(a-U),,,, 

bestimmt werden. 

Diese Aufgabe läßt sieb leicht mit Hilfe der Differenzen auf- 
lösen. 

Durch wiederholte Subtraktionen bestimmt man die Reihe der 
Werte 



A /(a), A f(a + Ä), A f(a + 2h), . . ., A f(a + n - Ih), 
A» f(a), A« /(ötH- ä), AV(a + 2h), . . ., AY(a + nT^h), 

A"-V(a), A-V(a + Ä), 

A- fia). 

Da aber A^fix) konstant ist (Nr. 4), so ist 

A"f(a + Ä) = A"/"(o). 
Mit Benutzrmg der Formel 

A*-V(« + Ä) = A*-V(«) + AV(aj) 
und durch wiederholte Additionen erhalt man die Werte 

A""" V(öt + 2Ä) = A"~ V(« + A) + A" f(a + h), 

AT'^fia + Sh) = A""V(ö^ + 2Ä) + AT'^fia + 2h), 



A* f(a + n- Ih) = A* f(a + n- 2h) + A» /(a + n- 2h), 

A f{a + nh) =A /(a + n-lÄ) -f A» ir(a -|- wT^l Ä), 

/■(a -f n -f lÄ) = /(a + wÄ) +A /*(a + wÄ). 

In derselben Weise werden auch die Werte 



fia + n + 2h), f(a + n + 3Ä), . . . 
berechnet. 
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um f(a — h) zu bestimmen, bemerken wir, daß 

A'fia -h)= AV(o) 
ist, und daß die Relation 

A*-V(a;) = A*~V(a; + h)- A''f(x) 
besteht. 

Wiederholte Subtraktionen geben uns die Werte 

A'-^fia - Ä) = A"-^f(a) - A" f(a- h), 
A^-^fia -h) = A"- V(o) - A"- V(ö - h), 

A« f(a - Ä) = A» f{a) - A» f(a - h), 

A f(a-h) = A f(a) - A« f(a - h), 

f(a -h)= . f(a) - A f(a- h). 

Ebenso findet man auch die Werte 

f(a - 2h), f{a - 3ä), . . . 
Als Beispiel nehmen wir 

a 1, Ä=l, f(x) = x^ + x*-2x-l. 

Die Rechnung gibt 

fi-l) = l, m 1, /-(l) 1. 

Durch Subtraktionen findet man 

Af{-1) 2, A/'(0) = 0, 

AY(- 1) = 2. - 
Nach der bewiesenen Formel (Nr. 4) ist aber 

A»f{x) = 6 
för jedes x. 

Die Anwendung der Additionen -gibt uns 

AY(0) = 8, A/(l) = 8, m=^l; 
A*/(l) = 14, A/(2) = 22, /(3) = 29. 

Um /■(— 2) und /"(— 3) zu bestimmen, sind folgende Rechnungen 
auszuführen: 

AY(-2) 4, A/-(-2) = 2, /•(-2) = -l; 

AY(-3) 10, A/-(-3) = 12, /•(-3) = -13. 



J « w - 



10 
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Diese Besoltate lassen sich in der Tabelle zusammenstellen: 



X 


/ 


A/- 


A»/ 


AY 


-3 


-13 


12 


10 


6 


2 


— 1 


2 


-4 


6 


-1 




-2 


2 


6 





— 1 





8 


6 


1 


9 


8 


14 


« r 


2 




22 






3 


29 









Aus den Vorzeichen von f(x) erkennt man, daß die Gleichung 
f(x^ = drei reelle Wurzeln hat, die in den Intervallen (— 2, — 1), 
(- 1, 0) und (1, 2) liegen. 

Diese Methode führt immer zur Absonderung der Wurzeln einer 
numerischen Gleichung, sobald der absolute Betrag der Differenz 
zweier Wurzeln kleiner als Eins ist. 

8. Belationen zwischen den sukzessiven Werten und 
Differenzen einer Funktion. Wir stellen uns die Aufgabe, die 
Werte 

durch 

f{d), f(a + h), /•(a + 2Ä), ..., f(a + nh) 

auszudrücken. 

Um die Schreibweise abzukürzen, wollen wir folgende Bezeich- 
nungen einfahren: 

f(a) = UQ, f(a + h)==Ui, f(a + 2h)=^u^, . . ., f(a + nh)=:Un. 

Die gesuchte Relation kann auf symbolischem Wege ge- 
fanden werden. 

In der Gleichung 

Awjfc = w*+^— w* 

können symbolisch die Indizes durch Exponenten ersetzt werden. 
Wenn man aus der Gleichung 

Am* = m*+^— w* 

den gemeinsamen Faktor u^ weglaßt, so wird 

A = M - 1. 



> < .c. .•- : 



!*''*• »,* • 
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Darans folgt 

Am ^m ^. ^ »» — 1 I w(w — 1) ^,m — 2 , / -t\m 

^u — mu H — . g u h (— 1) . 

Nach Multiplikation mit u* und Ersetzen der Exponenten darch 
Indizes erhält man 

Am . m(m — 1) I / •« \»» 

Uk=>UiJ^m—mUk^m^l+ \^ ' ^Jb+tn-2 1" (— 1) W*. 

Diese Formel ist damit zwar gefunden; aber noch nicht bewiesen. 
Das kann in verschiedener Weise geschehen. 
Wenn man in der Gleichung 

a durch a + h ersetzt^ so vermehren sich alle Indizes um Eins und 
es wird 

A Wjb-j-i = Uk+2 — w*+i. 

Durch Subtraktion findet man 

In derselben Weise fährt man weiter fort 

A*%+l = Wife4.8— 2Wä+2+ M*+l, 

A^Uk = (Mife+8 — 2Wife-|.2 + Wife+i) — (Uk+2 — 2w*_|.i + Uk) 
= UkJ^B — 3 Wife+2 + 3^*4-1 — Ujty 

A^Uk =(w*+4— 3%4.8 + 3wjfc+2— Wife+i) — (w*+s— 3wjfc4.2+3w*+i— wjt) 

Damit ist die Richtigkeit der gefundenen Formel für 

m = 1, 2, 3 und 4 

bewiesen. Durch vollständige Induktion läßt sich zeigen^ daß die 
Formel allgemein gültig ist. 

Wir werden aber einen andern Weg einschlagen. Aus der 
früheren Betrachtung folgt, daß A^Uk die Form hat 

weil wir immer mit linearen und homogenen Funktionen der u zu 
tun haben. Die Koeffizienten Aq, ä^, . . ., Am sind von der Fuuk- 
tioR f(x), von a und h unabhängig. Deswegen * genügt es, diese 
Koeffizienten in einem speziellen Falle zu bestimmen. 
Wir setzen 

a=rO, Ä=:l, f{x) = (f 
und finden 
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Folglich ist 

Die Koeffizienten J. sind aber ans dem binomischen Lehrsatze 
bekannt. 

Es ist also allgemein 

Setzt man in der bewiesenen Formel 

Jc = 0, m = 1, 2, 3, . . ., n, 
so erMlt man die gewünschten Formeln 

AMo=Wi-Wo; 

A*Mq = ^8 — 3t<j + 3^1 — Wo, 



An , W(n — 1) , / ^vn 

Wo=Wn— WWn-l+ '^ g ' Un-2 h (—1) % 

9. Jetzt wollen wir znr umgekehrten Aufgabe übergehen. Es 
sollen die Werte 

f(a), f(a + h), /•(a + 2Ä), ..., f{a + nh) 
durch 

f(a), Af(a), AY(a), ..., AV(«) 

ausgedrückt werden. 
Aus der Formel 

folgt 

Symbolisch ist 

oder 



Aw*=WA-f.l— w* 

M = 1 + A. 



Daraus ergibt sich 

w- =l + mA + ^?^^^A«+...+A", 

1 * A 

1 A I W (wt — 1) A 9 I I A m 

Wm+*=w*+mAw*+-Yr2-^A*M*H hA w*. 
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Um diese Formel zu beweisen^ begimien wir mit den Fällen 
w = 1, 2, 3. 

Wir haben schon gesehen^ daß 

Daraus folgt (Nr. 3) 
Awifc+i= Ama+ A*w*, 

w*+2= (w* + AWjfc) + {Auk+ A^Uk) = Wife+ 2AUk + A^w*, 

Aw*+2= Awa+2A*w*+ A^M*, 
Wa4.8= {Uk+ 2Aw*+ A^Uk) + (Am*+ 2A^Wi+ A^Mife) 
= M*+ 3Aw*+ 3A*Wjfe + A'wi. 

Wir haben hier immer lineare und homogene Funktionen von 
Uk, AUk, A^Uk, . . ., folglich ist 

Uk+m= B^Uk + JBi A w* + JBa A^Uk + h BmA'^Uk, 

wo die Koeffizienten von der Wahl der Funktion f(x) unabhängig sind. 
Für 

findet man 

c*+m= JBoC* + JBiC*(c - 1) + BjC*(c - 1)* + . • . + JBn,o*(c - 1)"» 

Setzt man c = 1 + < ein, so wird 

{l + tY=B^+ B^t + B^t^+ ^'^ + Bnt\ 
Die gesuchten Koeffizienten haben also die Werte 

Es bestehen daher die Formeln 

W3 = % + 2 A Wo + A * Wo 7 

«*8 = ^0+ 3Awo+ 3A»Wo+ ^'%; 



Wn = Wo + w ^ ^0 + ^1.2 ^* ^0 H f- A" Wo- 

Alle diese Formeln können durch die Formel 

U,= Uq+zAUq+ \^^ A^Uq+ h ^ 1.2-8...n ^ ^0 

umfaßt werden. Hier kann nur die Werte 

«™al,mor, XT = 0, 1, 2, 3, . . ., W 

annenmen. 

Damit schließen wir dieses Kapitel und gehen zu Anwendungen 
der bewiesenen Sätze über. 
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Zweites Kapitel. 
Interpolation. 

10. Ctonaue Interpolation. Die einfachste Aufgabe ist die 
folgende: Es sind die Werte einer ganzen Funktion f(() von nicht 
höherem als w*®"^ Grade gegeben 

Es soll der Wert der Funktion für t = x berechnet werden. Diese 

Aufgabe war schon gelöst für den Fall^ daß — ^— gleich einer ganzen 
Zahl ist (Nr. 7). 

Die bekannten Werte können wir wieder mit 

%> ^1^ ^; • • V w» 

bezeichnen. Es ist also der Wert der Funktion u» gegeben für 

;gf = 0, 1, 2, . . ., w. 
Nach der bewiesenen Formel (Nr. 9) ist 

w^=Wo+^^Wo+-Y72-A^Wo+*- + i.2.3...n ^ ^o 

für diese Werte von z. 
Wenn wir setzen 

" — r' 

SO ist 

, (< — a) (< — a — Ä) . . . (< — a — n — 1 Ä) A « 
"^ l-2-8...n ^ **o 

eine ganze Funktion, bei der die Relation 

^(0 = m 

für w + 1 verschiedene Werte von i besteht. Diese Relation muß 
identisch erfüllt sein, weil die algebraische Gleichung F{C) — f{i) = 
von nicht höherem als n^^ Grade nicht mehr als n verschiedene 
Wurzeln besitzen kann. Wenn wir t = x setzen, so finden wir 

/•(o.) = /•(«) + ^ A /•(«) + ^Il^^l^l^ AY(«) + •• • 

+ - 1.2.8. ..n -^ ^(«)- 

Das ist die Newtonsche Interpolationsformel, der wir schon 
früher begegneten (Nr. 4). 
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11. Andre Aufgaben der Interpolation wollen wir nnr für ganze 
Fnnktionen von nicht höherem als dem 4*®^ Grade behandeln. 

Es sind Werte einer ganzen Funktion f(i) vom 4*®^ Grade ge- 
geben 

f{a), f(a + h), fia + 2h), f(a + Sh), f{a + 4.K). 

Man soll für t = x den Wert der Funktion und ihrer Ableitungen 
berechnen. 

Nach der Newtonschen Formel ist 

m) = m + ¥ • ^m + ^*"t^VJ~'^ ^v(«) 

+ 123^» ^ ' W+ 1.2.3.4Ä* ^ ^W 

für jeden Wert von t. 

Der Einfachheit wegen setzen wir t = a + liZy die Funktion /*(<) 
verwandelt sich dabei in F(ji) und wir finden 



oder 



.(.-l)(.-2)(.-3) . 

^ 1 • 2 • 8 • 4 '^ IV*J 



+ (^-6.»+ll.»-6.).^. 
Durch Differentiation erhält man die Formeln 
F' (^) = A/-(a) + (2^-l).^ + (3.*-6^ + 2).^ 

i?"' (;?) = AY(a) + (g- 1) AY(a) + (6«*- 18« + 11). ^^^, 

F'" (0) = AY(a) + (2« - 3) ^^, 

F""(g)^A*f(a). 

Wenn e ^y für t = x ist, so sind die gesachten Werte 

f'(x)==iF'(f,), f"(x)^l,F"(y), f" (x) = ^^"' (yy 

rix)^^,F^(i,). 

Hier ist y durch — ^— zu ersetzen. 

Es seien z. B. für eine ganze Funktion f(t) von nicht höherem 
als 3*^ Grade die Werte 
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m^l, A2)-3, /•(4) = 53, /•(6) = 199 
gegeben. 

Es soll /""(l) bereehnet werden. 
In diesem Falle ist 

= 0, Ä = 2, x = l, y = Y 

r(l) = x[AY(0)-yAY(0)] 
A Y(0) = 53 - 2 • 3 + 1 = 48 , 
AY(0) = 199- 3.53 + 3-3-1 = 48, 

r(l) = ^(48-24) = 6. 

12. Für eine ga«iize Funktion 4*®"^ Grades und für h = l seien 
dip tV ftiH*e 

/•(O), A/-(0), AY(0), AY(0), AY(0) 
gegeben. Es sollen die Werte 

berechnet werden. 

Wenn wir die Bezeichnungen 

dfia) = f[a + ^) - /•(«), <JY(«) = Sf(a + ^) - df(a), . . . 
einführen, so ist nach der Newtonschen Formel 

/•g) = f(p) + mdfiO) + ?^^dY(0) + ^(^-^)(^-^) s'f(p) 

+ : — 1.2-3.4 ** ' W- 

Die Werte df(p), d'f(0), S^f(0) und S*f(p) lassen sich aas der 
Formel 

+ "'-;"r."r" ^v(o) 

bestimmen, die für jedes ^ gültig ist. Man findet 

KÄ) = m +ro^m-^ ^V(0) + ^ A« AO) - S ^V(0), 

^(Ä) = ^(0)+Ä^^(Ö)-ife^V(0) + jJ^AY(0)-iSAY(0), 

/•(i^) = A0)+Äa^(0)-|Lay(0)+J^a3/-(0)-S^V(0), 

^© = AO) + ^ A AO) - ^ AY(0) + ^ AY(0) - ^ AY(0). 
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Die Reihe von Subtraktionen gibt 

Sm = ^ A/-(0) - A AY(0) + ^ AY(0) - ^ AY(0), 

*/•© -ro^m-iö ^V(0) + äSö ^V(0) - S AY(o), 

*^(Ä) = lö ^ ^0) - 4 AY(0) + ^ AY(0) - ^ AY(0), 

«^/^O = ro ^f(^) - Wo ^V(0) + ^ AY(0) - ^ AY(0). 

*V(0) = iJö^V(0) - i^ AY(0) + iS AY(0), 
*V© = iJö ^V(0) - j^ AY(0) + ^ AY(0), 
*VS) = i^ ^ V(0) - ^ A»r(0) + ^ AY(0). 



*W) = j^AV(0)-^AY(0), 
*V© = i4AY(0)-^AY(0). 

Die gesuchten Formeln sind also 

*/(0) = ^ A/-(0) - 4 AY(0) + ^ AY(0) - S AY(0), 

*^V(0) = i^ AY(0) - i4 AY(0) + S^V(0), 
'^V(0) = i^AY(0)-^AY(0), 

Wir haben schon gesehen (Nr. 7), daß die Gleichui^ 
im Intervall (— 2, — 1) eine Wurzel besitzt. Die Gleichung 

fix) = g? (- 2 + a;) = 

hat aber eine Wurzel zwischen und 1. 

Seliwanoff, DifTerenzezireclmung. 2 



18 Erster Teil. Differenzen. 

Nach der angegebenen Methode können die Werte 

berechnet werden. 

Wir wollen voraussetzen, daß f[^^ und f[ ^^ ) verschiedene 
Vorzeichen hab^. Dann hat die Gleichung 9>(^) = eine Wurzel 
im IntervaU (- 2 + ^^ - 2 + ^^). 

Jetzt ist in 9>(^) die Substitution 

anzuwenden. Wir erhalten eine neue Funktion 

die aucli eine Wurzel zwischen und 1 besitzt. 

Nach demselben Verfahren findet man zwei Zahlen ^ tind ^^, 
für die /i(^) iind A( lo ) ^«™ch^®<i®''i® Vorzeichen haben. Die 
Wurzel der Gleichung <p{t) = Uegt also im IntervaU 

(-^ + ro + w' -2^^iö + ^)• 

Auf diese Weise kann man die Wurzel der Gleichung mit beliebiger 
Genauigkeit berechnen. 

13. Wir wollen jetzt folgende Aufgabe der Interpolation losen. 
Es ist eine ganze Funktion 4*®"^ Grades zu bilden, die für 

fia), f'ia), /■(« + Ä), /•'(« + h), f(a + 2h) 

vorgeschriebene Werte besitzt. 
Die Funktion 

hat die Eigenschaften 

F{a) = f(a), F{a + h) = f{a + ä), F(a + 2A) = /*(« + 2ä). ' 

Dieselben Eigenschaften bestehen auch für 

0{x) = F(x) + (a: — a) (a? — a — Ä) (a; — a — 2h)q)(x), 

wie auch die Funktion (p(x) gewählt wird. 
Aus den Forderungen 

*' (a) = f (a), «' (o + A) = f (« + h) 
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lassen sich die Werte 9?(a) und q)(a + h) berechnen. Deswegen ist 
die Funktion 

YoUsiÄndig bestimmt. 

Es sollen z. B. die Bedingungen 

/(0) = 1, /'(0) = 1, /•(!) = 2, /•'(!) = 3, /•(2) = 7 
befriedigt werden. In diesem Falle ist 

A/-(0)==1, A^(l) = 5, AY(0) = 4, 

F(x) = l+x + ^(^=il.4 = 2x'-x + 1, 

0(x) = 2!b'— X + 1 + x{x - l)(jc — 2)g>(x) 

= 2a^—x + l+ (3^— 3a^ + 2x)q>(x), 

^'(x) = 4a; - 1 + (a:» - 3«» + 2x)q)\x) + (3«» -6x + 2)ip{x), 

$'(0) = - 1 + 29(0) = 1, 9(0) = 1, 

*'(1) = 3 - 9(1) = 3, 9(1) = 0, 

9(a;) = l + a;-(-l) a; + 1, 

f(x) = 2a;*- a; + 1 + (a;«- 3a;*+ 2a;) (- a; + 1) 

= — a;* + 4a;» — 3a;* + a; + 1. 

Wir wollen nocli ein Beispiel behaadeln. Es soU eine ganze 
Funktion 3*^ Grades die Bedingungen erfüllen 

m=h /•(i) = i, /•'(!) = !, r(i)=4. 

Nach dem Taylorsclien Satze liat die Funktion 

F(x) = 1 + (a; - 1) + 4^^=^ = 2a;*- 3a; + 2 

die Eigenschaften 

2^(1) = 1, J?"(l) = l, ^"(1) = 4. 

Dasselbe gut auch für 

0(x) = 2a;*- 3a; + 2 + (a; - l)»9(a;) 

bei willkürlicher Wahl von 9(a;). Da aber nur eine Bedingung zu 
erfüllen ist, so können wir 9(a;) konstant voraussetzen. Es soll 
*(0) = 1 sein, folgüch ist 

2-9(0) = l, 9(0) = 1 
'^^ /•(a;) = 2a;*-3a; + 2 + (a;-l)«=a;»-a;*+l. 

Dieselbe Aufgabe -könnte man in andrer Weise lösen. 

2* 
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Den Bedingungen f(p) = 1, /*(!) = 1 genügt die Funktion 

0(x) = 1 + ic(a; — ^)<p{x) =■ 1 + (x' — !>o)g>{x). 
Die Ableitungen dieser Funktion 

*'(») = {a^-x)<p' (x) + (ßx - l)<p{x), 
9"(x)=(x'- x)g>"(x) + 2{2x - l)(p'ix) + 2(p{x) 
haben für a; = 1 die Werte 

*'(1) = 9(1) = 1, ^"(1) = 2y'(l) + 29(1) = 4. 

Daraus folgt 

9(1) = 1, 9'(1) = 1, 

(p{x) == 1 + (o? — 1) = aj, 

f{x) = 1 + (aj* — x)x = (x^—- x^+ 1. 

Gbnz ähnlicli kann man ganze Funktionen höheren Grades bilden^ 
die vorgeschriebene Werte 

f{a + h), f{a + h), r(a + Ä), . . . , /•(«^-^)(a + h), 

f(a + mh), f(a + mh), f{a + mJi)^ . . ., /^«m-i)(a + mh) 
besitzen. Der Grad dieser Funktion kann nicht die Zahl 

a + a^+ a^-\ |-a,„— 1 

übersteigen. (Vgl. Markoff, Differenzenrechnung, S. 1 — 6.) 

In den betrachteten Fallen der Interpolation wurde vorausgesetzt, 
daß die Werte der unabhängigen Yariabeln in äquidistanten Intervallen 
auftreten, weil nur in diesem Falle die Differenzen anzuwenden sind. 

14. Angen&herte Interpolation. Für Werte von t, die in 
einem Intervall liegen, wollen wir die gegebene Funktion f(t) an- 
genähert durch die ganze Funktion F{t) ausdrücken, die gewissen 
Forderungen genügt. Der Fehler dieser Annäherung soU abgeschätzt 
werden. 

Z. B. soll die ganze Funktion F{t) vom nicht höheren als 
w*«"^ Grade sein und die Gleichung 

erfttllen für 

< = a, a + Ä, a + 2ä, . . . , a + nh. 

Es ist uns schon bekannt, wie man eine solche ganze Funktion 
büdet (Nr. 10). 
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Es sei X eine Zahl; die im gegebenen Intervalle liegt. Wir 
werden näherungsweise 

setzen und es ist der Fehler dieser Annäherung abzuschätzen. 
Die Funktion 

fW - FW - Ä i.2.8...(n + l) 

bei beliebigem Werte von K verschwindet für 

t = a, a + h, a + 2h, . . . , a + nh. 

Wir wollen K so wählen, daß diese Funktion noch für t ^x 
verschwindet. Hier ist x eine Zahl, die von a, a -f Ä, a + 2ä, ..., 
a + nh verschieden ist. 

Die Gleichung 

m-m -Fit) - K ^'-''^^\-:-']:inl7~ — = Q 

hat dann mindestens n + 2 Wurzeln. Nach dem Bolleschen Satz 
hat die (w + 1)*® Ableitung von 9(t) mindestens eine Wurzel, die 
wir mit | bezeichnen wollen. Diese Zahl ist eine mittlere zwischen 
der größten und der kleinsten unter den Zahlen a, a + fih und x. 
Es ist aber 

folglich ist 

Es besteht also die Formel 

, (x—a){x—a—h). . .(a?--a— n— 1^) y^n^/ \.{x-a){x-a-h), . .(x-a-nh) i7^_|_i)/fc\ 
^ 1.2.8...n.Ä« ^ 'W+ 1.2.8...(n + l) f ^«''- 

Das ist die Newton sehe Interpolationsformel mit dem Best- 
gliede von A. Cauchy (vgl. Cauchy, Oeuvres, 1 serie, tome 5, 
p. 422 oder Comptes rendus, tome 11 [1840], p. 787). 

Besonders sind die Fällg zu bemerken 

„ = 0, fix) = f(a) + (x - a)^), 

n = 1, fix) = fia) + ^ Afia) + (^-«)(»-«-^) /•''(!), 

n = 2, fix) = fia) + ^ Afia) + (^-^X^-«-^) aY(„) 



+ 1T2T3 / ^«> 

Die erste dieser Formeln ist die Lagrangesche. 
Durch Weglassung der Restglieder entstehen Näherungsformeln, 
deren Fehler durch die Restglieder ausgedrückt sind. 



1 
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16. Wir wollen einen besonderen Fall der Interpolation betrachten. 
Für t ^^ X soll die gegebene Funktion f(t) angenähert durch eine 
solche ganze Funktion 3^^ Gfrades F{t) dargestellt werden, daß die 
Bedingungen 

J?(- 1) = /^(- 1), F{0)^f(pi r(p)^r(Ol F{l)^f{l) 

erfüllt seien. ^ 

Die Funktion F(t) hat die Form 

Fit) = fi-l) + t±^Aa-l) + ^^A^fi-l) + ^^±^l^Ä, 
wo die Konstante Ä sich aus der Bedingung 

bestimmen ^ßt. Wir finden 

I' (t) = Af{- 1) + I (2t + 1) AY(- 1) + f (3<» - 1)A, 

/•'(O) = Afi- 1) + 4 AY(- 1) - ^Ä, 

Ä = 6A/-(- 1) + 3AY(- 1) - 6/"(0), 

A = 6[/-(0) - /•(- 1)] + 3[/-(l) - 2m + /•(- 1)] - 6/-'(0). 

Nach Reduktionen ergibt sich 

A = 3/-(l) - 3/-(- 1) - 6/-'(0). 

um die Interpolationsformel mit dem Bestgliede zu finden, 
betrachten wir die Funktion 

^it) = m - lit) - K^^^^. 

wo die Eonstante K so gewählt ist, daß 9{x) gleich Null ist. 
Wenn wir die vier Zahlen 

— 1, 0, 1, ^x 
der Große nach ordnen, so entstehen solche vier Zahlen 

«i; ötj, aj, «4, 
die den Bedingungen genügen 

a^ < o^ < ttj < a^. 

Diese Zahlen genügen der Gleichung 9{t) = 0. Nach dem 
ßolleschen Satze hat die Gleichung 9^(f) =^ mindestens drei solche 
Wurzeln 6^, \f 63, daß 

% < 61 < »2 < 63 < ttg < 63 < a^. 
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Keine der Zahlen i^, b^ und h^ ist also gleich NulL Außerdem 
ist diese Gleichung noch für < = befriedigt. Die Gleichung *'(<) = 
hat also mindestens vier verschiedene Wurzeln. Nach dem Roll eschen 
Satze hat die dritte Ableitung von *'(<) oder 0^(t) mindestens eine 
Wurzel < = I, die zwischen a^ und a^ liegt. Da aber 

ist, so ist ^^\^ = r(t)-K 

Wir haben also die folgende Formel bewiesen 

/•(^) =. f(- 1) + ^ Are- 1) + <^±^ AY(- 1) + ^ ^"+/.^,".^r'^ 

"^ 12.34 ' ^^}' 

Der Koeffizient A hat den angegebenen Wert. 

16. In andern Fällen werden die Bestglieder mit Hilfe der- 
selben Betrachtungen gebildet. 

Wenn z. B. die Näherungsformel 

genau ist für ' ^ ^ ^^ 

t = aj a + A, a + 2A, a + 3A, 

f(t)^F\t) fär < = a, a + A, a + 2A, 

fn(t)^F"(t) für ^ = a, 

dann ist i^(0 eine ganze Funktion vom 7**^ Grade und die Näherungs- 
formel f{x) = F(x) enthält den Fehler 

(x — a)' (a? — a — Ä)* (a? — a — 2Ä)* (aj — a — 3 h) 



1-2.3-4Ö-6-7-8 



/(«)(!). 



Die Zahl | ist zwischen der größten und der kleinsten unter 
den Zahlen a, a + SA, x enthalten. (Vgl. Mark off, Differenzen- 
rechnung; S. 6—8.) 

17. Anwendung der* angen&herten Interpolation anf die 
Berechnung der Wnrzeln numerischer' Gleichungen. Wir 

machen die Voraussetzung, daß die numerische Gleichung f({) = 
im Intervall (a, b) nur eine Wurzel besitzt und daß /*"(<) immer 
dasselbe Vorzeichen in diesem Intervall behält. 

Da die Vorzeichen von f(a) und f(b) verschieden sind, so ist 

f(a)-fil>)<0. 

Wir werden f(t) naherungsweise durch eine solche lineare 
Funktion F(t) ausdrücken, ftir die 

F(a) = f(a) und F(b) = f{b). 
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Hier ist h = h — a und deswegen ist 

F(t) = fia) + ^ [/•(&) - /•(«)]. 
Für t = x wird 

Den Wert x wählen wir so, daß F(x) = wird. Folglich ist 

Es liegt also a; zwischen a und &. Für diesen x kann nicht 
/*"(!) = sein. Die Gleichung 

wäre dann befriedigt für t == a, x und b. Nach dem Bolleschen 
Satze hätte die Gleichung ^(f) = mindestens zwei Wurzeln zwischen 
a und b, und die Funktion 

müßte in diesem laterrall ihr Yorzeichen wechseln; was gegen die 
Yoraassetzung ist. 

Weil F(x) = ist, so ist 

Der betrachtete Wert von x läßt sich einfach durch Deter- 
minanten ausdrücken. Es ist 



X = 



m m 

a h 



m m 
1 1 



Wenn t im Intervall (a, b) liegt und 

fia)f\t)>0 

ist, 80 ist f{a)f"{l)>0 und foJ^Uch 

f{a)f{x)<0; 

da {x — a){x — b) negativ ist. Die Wurzel der Gleichung f(t) = 
liegt also zwischen a und x. 

Wenn wir a und 6 als erste Annäherungen der Wurzel be- 
trachten, so sind 

_ bm-am 

zweite Amüherungen. 

Es sei 

f(a)f"{t)<0, oder mf"(t)>0, 

für Werte von t, die im Intervall (a, b) liegen. 
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Dann ist 

/•(6)r(l)>0 und f(h)f{x)<0. 

Die Wurzel von f{t) = liegt zwischen b nnd x. Die zweiten 
Annäherungen sind in diesem Falle 

bm - am 



a, = 



~M^m ' ^^ - ^• 



In derselben Weise kann man von a^ und b^ zu den folgenden 
Äniiäherungen a^ und b^ übergehen usw. (Vgl. Weber, Lehrbuch 
der Algebra, I, § 116.) 



18. Anwendung der angenäherten Interpolation auf 
Berechnnng der Logarithmen nnd Antilogarithmen. Wir 

nehmen z. B. eine Tabelle fünfstelliger gewöhnlicher Logarithmen 
aller ganzen Zahlen von 1000 bis einschließlich 9999. Mit EUlfe 
der Werte log N und log {N + 1) soll log (N + x) berechnet werden, 
wobei < a; < 1 ist. 

Die Literpolationsformel 

■f{x) = m + X Am + "-^ra) 

gibt die Näheruugsformel ^ 

log {N+x) = log N+ X [log (N+H)- log JT] 

mit dem Fehler 

x(l — x) log« 



s = 



1-2 (JV+ D« 

Es ist 

e<3<yiÖ, loge<^, 

iV+|>10». 
Daraas folgt 

^<i6^« •'^«^ «<7.10-«. 
Dieser Fehler hat keinen Einfluß auf die fünfstelligen Logarithmen« 

19. Die Berechnung der Antilogarithmen führt auf die Frage: 
Es sind die Werte gegeben: 

logJV^-a, log (JV^ + 1) - 6. 

Aus der Gleichung • 

log (N+y)^a + x 

soll y berechnet werden für einen Wert von x, der zwischen und 
b — a liegt. 



V 
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Wir setzen ^, . ^^ ^ ^r. i 

f(x)=-N + y-=^ 10» + * 

Die Anwendung der Interpolationsformel 

fix) = f{0) + ^ A m + ?(^ZM£) fn(^^ 
gibt die Näherung 

mit dem negativen Fehler 

s = ~ ^^^"'^^"^^ »10^ + €(1 10)1* 
Es ist 

Aus der Relation 

folgt ^""=10 

110 = ,— - 
löge 

Da aber 10<' + «<10» ist, so ist 

Wenn wir bemerken, daß 

10» = iV+ 1, & - a = log (N+ 1) - log 2f = log(l + ^), 
so finden wir 

-«<|(^-+ l)-aoge)-»[log(l +^)J- 

Dieser Ausdruck kann noch vereinfacht werden. 
Für positive ist 

(p(z) = log (1 + ;8f) — i8? log e < 0, 

weil 9(0) = und 

sind. Deswegen ist 

log(l+;^)<|^loge 
und folglich N4-i 

— F ^ — — • 

Da aber N+l£10^, N'^ 10^ so ist die positive Zahl 

^ -«<2.10-8. 

In diesem Kapitel ist die Interpolation nur als Anwendung der 
Differenzenrechnung behandelt worden. Ausfuhrliche Mitteilungen 
über diesen Gegenstand findet man im Referat von J. Bauschinger 
(EncykL der Math. Wiss., I, S. 799—820). 

* 1 10 ist der natürliche Logarithmns von 10. 
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Drittes Kapitel. 
Angenäherte Berechnimg bestimmter Integrale. 

20. Einleitende Bemerknngen. um das Integral 1 f(x)dxy 



a 



WO a<.b ist; angenähert zu berechnen, gebraucht man eine der 
Interpolationsformeln (Nr. 14 — 16) 

fix) = F{x) + <p(a;)/'(i»(|). 

Hier sind F(x) und g)(x) rationale ganze Funktionen und | ein 
Mittelwert. 

Nach Integration findet man die Näherungsformel 

b b 

Jf(x)dx==jF{x)dx 

a a 

mit dem Fehler 



o 

^fg>{x)f(p){i)dx. 



a 



Es ist I eine unbekannte Funktion von X] wir wissen nur^ daß 
für a<x<ib auch a < | < 6 ist. Der Fehler der Näherungsformel 
ist nur dann abzuschätzen, wenn die Funktion (p(x) für alle x zwischen 
a und b ihr Vorzeichen behält. Nach einem Satze der Integral- 
rechnung ist dann 



= r''Kn)J'Pix)dx, 



a 



WO a<iri<b ist. 

Die Aufgabe ist also auf die Integration ganzer rationaler 

Funktionen zurückgeführt. Die Bechnungen sind noch einfacher, 

wenn die Grenzen des bestimmten Integrals — 1 und + 1 sind, z. B. 
1 



— 1 



GPo + Ä^ + i>2^^ + i>8^ + PA,^)dX = 2|?o + TjPs + TJP4- 



Um die Grenzen des gegebenen Integrals in — 1 und + 1 zii 
verwandehi, setzen wir 

a? = ^(0 

und es sollen die Bedingungen erfüllt werden 

^(~l) = a, ^(1) = 6. 
Nach der Newtonschen Interpolationsformel ist 

I *+l /x \ h — a. , a + 6 
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Vermöge dieser Substitution findet man 

b 1 



Jf{x)dx ^ ^fHi)dt 



—1 

In den folgenden Paragraphen werden wir verschiedene Inter- 
polationsformeln anwenden. Daraus entstehen verschiedene Methoden 
zur Berechnung des gegebenen Integrals. 

21. Methode der Bechtecke. In der Formel 

b 1 



Jf{x)dx==^-^p{i)dt, 



— 1 

setzen wir (Nr. 14) 

•f (0 = J=(- 1) + (^ + 1) ^'{ty), (-l<ty<t). 

Daraus entsteht die Näherongsformel 



/ 



1 

^(t)dt == 2/(- 1) 



mit dem Fehler 

1 



—1 _1 

Hier ist — 1 < r < 1. 

Um die Funktion f(x) einzuführen, bemerken wir, daß 

•f'(0=/'(a')-^- 

Wir finden die Näherungsformel 

b 



Jf(x)dx = (6 - a)f{a) 



mit dem Fehler 

wo a < I < 6 ist. 

Geometrisch ist das gegebene Integral eine Fläche, die in der 
Näherungsformel durch ein Eechteck ersetzt wird. Daher rührt der 
Name der Methode. 

Um den Fehler zu vermindern, zerlegen wir das Intervall (a, b) 
in n gleiche Intervalle. Das Integral verwandelt sich in die Summe 

b a-\ a+2 b 

f(x)dx ==lf(x)dx + I f(x)dx + '" + lf{x)dx, 

a a a-\ a+w — 1 

n ' n 
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Nachdem wir auf jeden Snmmanden die Nähemngsf ormel angewandt 
haben; finden wir 



/' 



f{x)dx = ^ (yo + ^1 + y« + • • • + y»-i)> 



a 

WO der Kürze wegen die Bezeichnung 

eingeführt ist. 

Der Fehler dieser Näherungsformel ist 



6 = 



1 (6~a)\.., 



2 n 



^[/•'(ii)+/"(y + •••+/•'(!«)], 



wo a + (ÄJ — 1) < It < a + Ä ist. 

* 

Diesen Ausdruck kann man vereinfachen. 

Es wird vorausgesetzt, daß f{x) endlich und stetig im Intervall 
(a, b) ist. Diese Funktion hat dann eine obere Ghrenze Ä und eine 
untere Grenze Ä* 

Aus der Reihe von Ungleichheiten 

^ > /"'(ii) ^ J5, 4 ^ /•'(!,) ^B,...,Ä^ f'iL) ^ B 
folgt 

nA ^ /"(lO + /•'(!,) + ... + /•'(!„) ^ «5 

"^^ ^ > f(gi)+/"(g,)+--+r(£«) > _ß 

= ^ = • 

Nach einem Satze von Weierstraß kann f(x) jeden Wert 
zwischen A und B annehmen. Deswegen gibt es eine solche Zahl |, 

^^^"^ r(Si)+r(s,)+-"+r(§n) _ ^t ^j^y 

Es ist also 

Wenn wir die Interpolationsformel 

Sit) = /(i) + (< - i)j:'(o 

benutzen und. dieselben Rechnungen wiederholen, so finden wir die 
Näherungsformel 



/' 



f(x)dx = ^ (yi + y, + • • • + 2^«) 



a 

mit dem Fehler ^ ,, ,, 



^x=^i^^^Ai)- 



^Tannery, Introdnction ä. la th^orie des fonctions d'une variable, 
p. 101 et 126. 
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Wenn f{po) immer dasselbe Vorzeichen im Intervall (a, 6) hat, 

h 

Jf{x)d: 



SO ist der Wert von / f(x)dx zwischen 



a 

6 — a 



(yo + yi + ys + '-' + j/n-i) 



n 
und 

-^ (y 1 + ^2 + % + ••• + y«) 

enthalten, weil s und a^ verschiedene Vorzeichen haben. 
22. Um das Integral 

6 1 

Jf{x) dx = ^ A (0 ^^ 

a — 1 

ZU berechnen, kann man nicht die Formel 

/(O = /(O) + tP{t^) 

gebrauchen, weil t, der Paktor bei J^^it^y das Vorzeichen im Intervall 
(— 1, 1) wechselt und der Fehler der Näherungsformel nicht ab- 
geschätzt werden kann. Es ist zu setzen („Markoff^^, S. 53) 

und man findet 

1 

/(Ot?^ = 2/(0) + ^J="(r). 
—1 

Daraus folgt die Näherungsformel 



ß 



ff(x) -(h- a)f(i±^) 



mit dem Fehler 

2i(&-«)r(6)-2 YTäTs TT2* 

um eine genauere Formel zu bekommen, zerlegt man das Intervall 
(a, &) in n gleiche Intervalle. Wenn man dieselbe Methode wie in 
Nr. 21 benutzt, so findet man die Nähemngsformel . 

h 



ß 



f{x)dx — (yi + ys + • • • + yan-i) 



SS s 

a 

mit dem Fehler ^ ^ ^ (6_a)» ^ fu^ 

* ~ 2n» * 1 • 2 • 3 * 1 • 2 * 
Hier ist wieder 
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28. Methode der Trapeze. Znr Bereclmung des Integrals 

b 1 

a —1 

werden wir die Funktion /(^) angenähert durch eine lineare Funktion 
9(ß) ausdrücken^ die den Bedingungen genügt 

*(-!) = /(-!), «(1) = J=(1). 
Dadurch entsteht die Interpolationsformel 

j:(<) = /(- 1) + *-±^ \f{l) - H- 1)] + (*+;)(^^-^) j:»(^). 

Die Anwendung derselben gibt 

1 

//(Od< = /(-l) + /(l)-|^"(r). 

— 1 

Daraus folgt die Näherungsformel 

b 



Jf{x)dx = *^ \f{a) + fm 



mit dem Fehler 

1.23 * 1-2* 

Die rechte Seite der Näherungsformel ist der F^cheninhalt 
eines Trapezes. 

Nach Zerlegung des Integrals in n Summanden findet man die 
Näherungsformel 

b 
a 

oder 

b 

Jmdx = ^ (yo + 2yi + 2y, + . . • + 2y,_, + y„) 

a 

mit dem Fehler 

_ _!_ (6-a)» r(S) ^ 
^ n* 1.2-2 12 

Wenn f\x) dasselbe Vorzeichen hat für alle x zwischen a und 
6; so haben s^ und b (in Nr. 22) verschiedene Vorzeichen. Deswegen 



ist das Integral f f(x)dx zwischen 
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h — a 



n 



(yi + 2/3 + 2^5 H hVin-i) 



2 S 2 2 

und ^_^ 

-2^ (yo + 2^1 + 2ya + • . • + 2y„ _ j + y„) 
enthalten. 

24. Slmpsonsche Formel. Wenn wir die Interpolations- 
formel (Nr. 15) 

S(t) = .f (- 1) + ^' A /(- 1) + ^* A*/(- 1) 

. (t + i)t(t-i) {t + i)t'(t^i) i^, . 

+ ^ iT2T8 + 12.84 "^ W 

anwenden^ so wird 

\t)dt ^ 2/(~ 1) + 2A/(- 1) + 4 A«/(- 1) - i^ J^^{r). 



— 1 

Es ist aber 



A/(- 1) = m - /(- 1), AV(- 1) = /(l) - 2/(0) + /(- 1), 

^^(0 = rix) ■ (^)*. 

Nach einigen Reduktionen finden wir die Simpsonsche Formel 



/' 



««)* - ^ [/•(») + 4,'(i±i) + m] 



a 

mit dem Fehler (6-^)^ fiv(g) 



1-2-8.4-5 1.2-8-4 

Wenn f{x) eine ganze Funktion vom 3*«^ oder niedrigeren Grade 

ist; SO ist das Bestglied gleich Null und die Simpsonsche Formel 

ist ganz genau. Z. B.: 
1 

(2x + 3) (2a; - l)(2a: + b)dx == -1 [- 3-5 + 5 • 7] = y. 



Wenn aber f{x) eine andre Funktion ist, so kann man den 
Fehler verkleinem; nachdem man das gegebene Integral in n Sum- 
manden zerlegt hat. Man findet dann die Näherungsformel 

b 



ß 



a 

mit dem Fehler 



f(x)dx =^-^{y^ + 4j/£ + 2^1 + 4y^4- • • • + 4y 2n-i + y«) 

8 8 2 



. 1 (6-ar A^(S) 



«* 1.2-8-4.5 1-2-3.4 

In diesem Kapitel haben wir mit geringen Änderungen die 
Darstellung vom Herrn A. A. Markoff angegeben. (Difierenzen- 
rechnung, S. 50 — 60.) 

Wir schließen damit die Anwendungen der Sätze über Differenzen 
und gehen zu den Summen über. 



Zweiter Teil. 

S n m m e n. 



Erstes Kapitel. 

Unbestimmte und bestimmte Summen. 

26. Defiliitionen. unter Summe von f{x) versteht man eine 
Funktion q>{x)y deren Differenz gleich f{x) ist, also 

q>{x + Ä) ~ (p{x) = fix). 

Diese Gleichung bedeutet, daß der Ausdruck 9(0; + A) — (p{x) 
so transformiert werden kann, daß er in f{x) übergeht. Dieselbe 
Relation wird auch bezeichnet durch 

q>{x) =^ 2f{x). 

Es ist vorausgesetzt, daß f{x) nicht in der Form q>{X'\-}i) — q>{x) 
gegeben ist. Sonst wäre 

q){x) = 2[(p{x + Ä) — (fix)] 

eine reine Identität, die keinen Nutzen haben könnte. 

Wir werden untersuchen, ob es möglich ist, daß eine andre 
Funktion ^(rz;) dieselbe Eigenschaft 

il>{x + A) - il>{x) = fix) 

hat. Wenn wir il)ix) — (pix) mit Fix) bezeichnen, so ist 

il>ix) = (fix) + Fix) 
und nach Nr. 3 

A^(a;) =* Aq}ix) + AJ?(a;). 
Es ist also 

AFix) =. oder Fix + A) = Fix). 

Beide Funktionen ij^ix) und q)ix) unterscheiden sich also um 
eine periodische Funktion mit der Periode A. 
Ist ^ix) eine beliebige Funktion, für die 

cP(a; + A) = cP(ic) oder A*(a;) = 0, 
so ist 

A[(pix) + *(ic)] = Affix) + A 0ix) ^fix). 

Seliwanoff, Differenzenrechnimg. 3 
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Also haben alle Funktionen^ deren Differenz gleich f{x) ist^ die Form 

wo 9{x) eine willkürUche periodische Funktion mit der Periode Ä ist. 
um die Aufgabe zu vereinfachen, werden wir von nun an voraus- 
setzen; daß X nur Werte der Reihe 

a, a + h, a + 2h, . . . 

annimmt. Dann ist 9{x) eine willkürliche Konstante, die wir mit G 
bezeichnen werden. Es ist also 

2f{x) »= g>(x) + a 

Dieser Ausdruck heißt unbestimmte Summe, weil der Wert 
von G ganz beliebig ist. 

26. Bnmmation einflEtoher Funktionen. Aus jedem Satz 
über die Differenz (Nr. 2 — 6) können wir einen Satz über die Summe 
herleiten. 

Aus der Relation Ax = h folgt Ay- = 1 und 



Weü 



2:i = J + a 



Ax(x — h)(x — 2h) . . . (a? — n — lÄ) = nhx(x — ä) . . . (a: — » — 2h) 
ist, so ist 

A -^xipc — Ä) (ic — 2Ä) . . . (a; — w — lÄ) = x(x — A) . . . (a; — n — 2ä) 
und 
Sx{x'-h){x-2h) . . .{x-n — 2h) ^ ^x{x-h) . . .{x ■--n-lh) + G, 

Damit man unter dem Summenzeichen n Faktoren hat, ist es 
zweckmäßig, n durch n + 1 zu ersetzen. Man findet 

Zx{x — A) (a? — 2A) . . . (a; — w — lA) ^r-Trih^^^'^ A) . . . (a; — nh) + (7. 
Wir hatten die Formel (Nr. 2) 



Im Nenner der rechten Seite sind mindestens zwei Faktoren^ 
weil w ^ 1 ist. 

Aus dieser Formel folgt 

^ 1 _ 1 1 . ^ 

x{x-\-h).,.{x + nh) nh x{x-\-h),,,{x-\-n-lh) 

Die Funktion unter dem Summenzeichen hat im Nenner mindestens 
zwei Faktoren. 
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Wenn wir w — 1 statt n einsetzen; so entsteht die Formel 

2 — =- r^^T^ — =^ + a 

x(x+h)...(x-\-n-lh) (n-l)Ä x{x+h) . ,.(x+n-2h) 

Euer ist n > 1. Dem Wert n = 1 entspricht eine besondere 
transzendente Funktion, die wir mit 2— bezeichnen. 

' X 

Die Relation 

Am* = m^ivnh — 1) oder A —. — - == m* 

gibt die Formel 

Sm' = -^ + G. 

Es ist, wie wir wissen (Nr. 6), 

A COS w = — 2 sm -r- sm li* + -r-j- 
Wir setzen . ^ 

und wählen a und /3 so, daß 

, Au 



wird. Ans der Bedingongsgleichnng 

I 

2 






folgt Ä 

Deswegen ist 

A cos (a; — 2^) = — 2 sin — • sin x. 

Daraus folgt ^ 

co8(a:--) 



und 



2 sin — 
2 



cos 



= sma? 



(-4) 



2? sin a; = ^ h (7. 

2sm- 

Ganz analog findet man die Formel^ 

sin(a;--) 

2;C0SiC^ ^ I. + g. 

2sin — 
2 

27. Snmmation aii8ammenge8etzter Fimktioneu. Wenn 
A ein konstanter Faktor ist, so ist 

^A^^{x) = ÄAq)(x) 
^^'^'^ 2[ÄAg>(x)]^Ä(p{x) + a 
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Wenn wir A(p(x) mit f(x) bezeichnen^ so entsteht die Formel 

2Äf{x)^A2f(x), 

Die Konstante C ist schon im Ausdrucke 2f(x) enthalten. 
Wenn wir eine Reihe von Gleichungen haben 

so ist 

Daraus folgt 
^\fi{^) + U{x) +•'• + U{x)] ^ ip,{x) + (p,(x) + ' • ' + (p„,(x) + G 
oder 

2r/i(^) + f2{^) +'"+ u{^)\ - ^fi{^) + ^u{^) + • • • + ^/-(4 

Mit Hilfe dieser Formel ist es leicht^ ganze rationale Funktionen 
zu summieren. 
Es sei 

fw - A + A^- + A- — -JT2 +"' 



, j (a;--a)(rc--a--Ä)(rc--a — 2Ä) . . .(rc — a — n — lÄ) 
■^^ 1.2.3...n 

Dann ist 

■• '■^ 1.2.3...(n + l).Ä 

Die Summe einer ganzen Funktion n*^ Grades ist also 
eine ganze Funktion vom (w + 1)*®^ Grade. 
Für Ä = 1 ist z. B. 

a;^ ^ a? + a;(ic — 1); ic* == ic + 3a;(a; — 1) + a;(ii; — 1) (a? — 2) 

und folgUch 

2?a;« = C + {«(a; - 1) + \x{x - 1)(« - 2), 

Sa? = C + ^x{x - 1) + «(« - l)(a; - 2) + ja;(a;- l)(a;- 2)(a;- 3). 
Es ist aach leicht, gebrochene Funktionen von der Form 

m 

{x -{- a)(x -{- a -{-h) . . . (X'{-a-{-nh) 
zu summieren; wo' f(x) eine ganze Funktion ist. 

28. Partielle Summation. Aus der Formel xj^^ 

A[(p{x) • if{xy] = (p(x) • ^i;(x) + (p(x + h) • Aq>{x) ^^ 

^ ^ 2\ip{x) . t^^{x) + ^(a? + Ä) • t^q>{x)] = g>(a;) • ^(a;) + 
oder 
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Diese Gleichung liefert die Formel der partiellen Summation 

2J(p(x) • Af{x) = (p(x) ' ilf{x) — 2]f(x + h) • A(p{x). 

Die willkürliche Konstante erscheint, nachdem die Summation 
auf der rechten Seite ausgeführt ist. Mit Hilfe dieser Formel kann 
man das Produkt (p(x) • f{x) summieren, wenn man einen Wert von 
2^f(x) kennt, und wenn die Operation 2]ilf(x + h) - A(p(x) ausfahrbar 
ist, wo ^(a;) = 2Jf(x) ist. 

Wir wollen z. B. US^- x für ä = 1 ausrechnen. 

Da ein Wert von 2? 3* gleich — 3* ist, so ist 

2:3*. a: = 2;a; • A (^ - 3-^) =- ic . 1^ . 3*- 2;| • 3'+^ • Ax 

Es wäre nicht zweckmäßig, die Summation beim Faktor x anzufangen. 
Da Ux = ^ r ist, so wäre dann 

= |3*a;(a? - 1) - 2](x + l)x • 3% 

und die gegebene Summe würde auf eine kompliziertere zurückgeführt. 
Noch ein Beispiel: 

— cos Ix —'i) 



Ux sin o; = 2x • A 



-cos(a.-|) 



2 8m 



= X 



2 sin — 



2 J 2 



2sin- 2sm- ^28m-j 

28. BigenBchafben der bestimmten Summen. Wir haben 
schon gesehen (Nr. 25), daß, wenn 

(p(x + Ä) - (p(x) = f{x), 11^(X + Ä) - if(x) = f{x) 

und X eine Zahl der Reihe . 

a, a + h, a + 2ä, . . . 

ist, (p(x) und ilf{x) sich durch eine Eonstante unterscheiden 

(p(x) == if(x)+ G. 
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Subtrahiert man voneinander zwei Werte dieser Funktionen för 
x = a + nh und x=^a + mh, so föllt C weg. Man findet die Gleichung 

q)(a + nh) — (p(a + mh) == fia -\- nh) — ^(a + mh). 

Der Ausdruck (p(a + nh) — (p{a + mh) bleibt also ungeändert, 
wenn man (p(x) durch einen andern Wert von 2f(x) ersetzt. Dieser 
Ausdruck heißt bestimmte Summe und wird bezeichnet durch 

a -\-nh 

^f(x) = (p{a + wÄ) — (p{a + mh). 

a-{-tnh 

Die Zahlen a + mh und a + nh werden Grenzen dieser Summe 
genannt. 

Aus dieser Definition folgt' 

2m = 0, 2f{x) = f{a + mh). 

Wenn c eine Zahl a + Tch ist, die zwischen a und h ^ a + nh 
liegt, so ist 

b 

^fia>) = 9(6) - ip(a) = [9(6) - 9(0)] + [ip(e) - 9(0)] 
oder beb 

a a c 

Wenn man diesen Satz mehrmals anwendet, so findet man 

b a+h a + 2ft b 

-2f(.'') =2/"(*) +2/"(^) +-"+2/"(^)- 

Daraus folgt eine wichtige Formel 
f{a) + f(a + Ä) + /•(«+ 2ä; + ... + f(b - Ä) = (p{b) - (p(a).] 

Hier ist /"(o?) die gegebene Funktion, ,(p(x) irgend eine Funktion, 
deren Differenz gleich f(x) ist, -^— eine ganze positive Zahl. Diese 
Formel kann nur dann nützlich sein, wenn (p{x) nicht die Form hat 

f(a) + f(a + h) + f(a + 2h) + '" + f(x - h). 

30. Berechnung einiger bestimmten Sammen. 1) Es sei 

s == a + {a + h) + (a + 2h) H + {a + n — Ih). 

Wenn wir a + nh mit h bezeichnen, so wird 

b 
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Ein Wert der nnbestiininteiL Summe I]x ist (Nr. 26) 

deswegen ist 

oder nach einigen Redaktionen 



n 



s = -^[a + {a + n - Ih)]. 

' 13 



2) 1.2-3 + 2.3-4 + ..- + 10.11.12-2^(^- 1)(^~2) 



= 1.1312.1110. 

4 



28 
^ 1 . S "*■ 8 . 5 "^ 6 • 7 "' '■ 21 . 23 j^as(x + 2) 2 46 23* 



X '(' 



7» 



nA 



5) 1 + cos h + cos 2Ä H h cos (w — 1)ä =^ cos o? 



. 2« — 1- , . h 
sin — - — Ä + Bin- 

2sin-r- 
2 

Wenn wir von beiden Seiten — abziehen^ so finden wir 

. 2n-l, 
1 sin— ^— Ä 

2^ + cos Ä + cos 2Ä + • • + COS {n — 1)ä = j- — 

2sin — 
2 

«Ä Ä 2n — 1, 

J*"- cos — — cos — - — h 

6) sinÄ + sin2ÄH h sin (w — 1)^ =^,sina; = ^r 

2sin — 

2 

7) Wir wollen jetzt mit Anwendung von Formeln in Nr. 27 die 
Summen von Potenzen natürlicher Zahlen berechnen. 

X 



X 



x{X'-i)(2x-l) 
6 
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X 

0»+ 1»+ 2»+ ••• + (a: - 1)» = '^a?=^^x{x-l) + x{x-l){x-2) 



Daraas folgt die bemerkenswerte Relation 
1»+ 2»+ 3»+ ••• + («- 1)»= [1 + 2 + 3 + ••• + («-!)]».* 

Wenn wir :;— r und statt iZJ* und a? summieren^ so finden wir 



X 



^T^ ~~ 123 ¥'lT2"^12^ 1-2-8 ' 



-^1.2.3""l.2.34 2 'l.2.3"*"l2*1.2~ 1.2.34* 



Im nächsten Kapitel werden wir die Funktion betrachten 

-2j1.2.3...(«-1) ^ 9^'»(^)' 

diC; abgesehen yon einem konstanten Faktor^ von Jacob Bernoulli 
eingeführt ist. 

Zweites Kapitel. 
Die Jacob Bernonllische Funktion. 

31. BestJmmnng der Koeffizienten. Nach dem in Nr. 27 
bewiesenen Satze ist die Jacob Bernonllische Funktion 

eine ganze Funktion n*^ Grades^ die wir in der Form 

schreiben werden. Das konstante Glied fehlt^ weil 9>n(0) = ist. 
Denn in Nr. 29 haben wir gesehen^ daß die bestimmte Summe gleich 
Null ist^ wenn ihre Grenzen einander gleich sind. Wir wollen die 
Koeffizienten A^, Ä^, A^, . . . , An—i bestimmen. 

* Vgl. Encykl. der Math. Wißs., Bd. I, p. 679. 
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Da 

X-{-l . OB . g-f-l . 

q>n{x + 1) =^1.2... (n~l) ^^1.2... (w-1) +^ 1.2...(w-l) 
ist, so besteht die Belation 

(3) Vn{x + 1) -^ g>„(a?) = i.2'^^(^_i) > 

die nur für ganze positive x bewiesen ist. Die ganze Funktion 

F{X) ^ ip„(x + 1) - (pnix) - .^.a"^. (^_i) 

vom nicht höheren als w*®^ Grade verschwindet für jeden ganzen 
positiven Wert von x. Also ist F(x) identisch gleich Null und die 
Relation (3) besteht für jedes x. 

Wenn eine andre ganze Funktion ^(i») die Eigenschaften 

i>ix+l)-tix) = ^^^'^~^^_^y V(0) = 

hat, so ist xl;(x) mit q>n(x) identisch, weil beide Funktionen vom 
n^^ Grade (Nr. 4) sind und außerdem die Gleichungen 

*(0) = 9>(0), ^(1) = g> (1), ^(2) = q>{2), . . . , ^(n) = ip(n) 

bestehen. Die Werte q>(l), 9?(2), . . . , (p(n) bestimmen sich sukzessive 
aus der Gleichimg (3). 

Durch die Differentiation von (3) findet man 

fp'n(x + 1) - fp^nix) = i.2'^,(^-2) ' 

Nach der Gleichung (2) ist (pn(P) = Än—i. Deswegen ist 
(f/n(x) — Afi—i eine ganze Funktion, die für o? = verschwindet und 

deren Differenz gleich r— ^ 7~zr^ ^^*- Folglich ist für w > 1 

oder 

(4) (pn(x) = (pn-l(x) + An-1. 

Aus dieser Eelation schließen wir, daß 

y„-i(a;)==A i.2!^.(^_l) + A i.2!^.(n-2) + - • + -^"-'^ 

ist, und daß die Koeffizienten J^, A^, . . . , -4«— g hier dieselben sind, 
wie in der Gleichung (2). 

Wenn wir in (3) x=^0 setzen, so finden wir (fnQ-) = oder 

(5) — ^2 1 ^1 L_ :^ L u^"-^-lA.-i = 0. 

^^^ 1.2...n^ 1.2...(w-l)^ 1.2...(w-2)^ ^ 1-2 ^^^-^ 
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Diese Relation gibt die Möglichkeit^ ^— i zu berechnen^ sobald 
die Koeffizienten Aq^ A^, . , ,, ^n~s bekannt sind. 
Aus der Definition 

X 



folgt, daß ^. j 

ist. Deswegen ist 

-4^== 1, -4^= — — • 

Die Gleichung (5) für w = 3, 4 und 5 gibt 

12. 34 + iT2T3 + TT2 + ^8 =^ ^' 
l-2.3.4.5"^1.2.3.4"^l-2-3^1.2^* 

Daraus findet man 

In dieser Weise kann man beliebig viele Koeffizienten berechnen. 

32. Eigenschaften der Koeffizienten. Wenn man in der 
Relation 

%i+lO- + «) - Vst+i W = i.2.8...2ife 



««* 



X durch — X ersetzt, so findet man 

(6) y»+i(i - ^) - y»+i(~ ^) ° 1 • 2 • 8 ■ ■ ■ 2ifc • 

Wenn man die Bezeichnung 

- fik+iO^ - «) = -2^(«) 

einfuhrt, so ist 

- y^+iC- *) = - 9'"+iCi - («^ + 1)^ = ■^('*' + !)• 

Die Gleichung (6) nimmt die Form an 

Fix + 1) - F(x) ^ ^:^^. 

Da aber F(0) = ist, so ist nach dem in Nr. 31 bewiesenen 
Satze F(x) = 92*+i(^) ^^®^ 
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Vergleicht man die Ableitungen der beiden Seiten dieser Gleichung, 
80 findet man mit Berücksichtigung der Gleichung (4) 

oder 

(ß) 9^2*(l~^) = 9>2ife(4 

Für Ä = 1 ist die Ableitung yon 9?3t(a;) keine Jacob Bemoullische 
Funktion, weil die Definitionsgleichung (1) in Nr. 31 für n = 1 keinen 
Sinn hat. Aus der Differentiation von (8) folgt für Ä > 1 

- [9^2*-i(l - ^) + -^a*-!] = 9>2*_i(^) + ^2*-l 

oder nach (7) 

— -42* — 1 = -^2*— 1. 

Es ist also 
(9) Ä2jt-i = für Ä > 1. 

Wir haben schon in Nr. 31 gesehen, daß Ä^ nicht Null ist, und 
daß -4^ = ist. Jetzt überzeugen wir uns, daß -^5 = 0, -4^ = 0, . . . sind. 

33. Das Torzeichen von 9>nix) im Intervall (O, 1). 

Wenn man in (7) 00 = -^ setzt, so findet man 



fp»+i © = 0- 



Wir wollen untersuchen, ob für Ä> 1 die Gleichung V2k+i(^)'^^ 
noch eine Wurzel im Intervall (0, 1) hat. 

Es sei vorausgesetzt, daß diese Gleichung außer und 1 noch 
die Wurzeln a^ und a^ besitzt, die den Ungleichheiten 

< ai < »2 < 1 

genügen. Euer darf a^ oder a^ gleich — genommen werden. Nach 
dem Eolleschen Satze hat ^^kA-ii^) ^ ^ mindestens drei Wurzeln 
\^ &2 ^^^ ^3^ ^® ^^^ Ungleichheiten genügen 

• * • 

< &! < a^ < &2 < % < 63 < 1. 

Die Gleichung ^'gjfc-i-i^^) "^ ^ ^^* niindestens zwei Wurzeln c^ 
und c^, die zwischen und 1 liegen. 
Es ist aber 

weil -42 t— 1 = ist. Die Gleichung ^^j*— 1(^) ^ ^ ^^^ außerdem die 
Wurzeln und 1. 

Aus der Voraussetzung, daß die Gleichung ?P2*+i(^) ^ ^ zwischen 
und 1 zwei Wurzeln hat, folgt also, daß die Gleichung ?P2*-i(^)'^^ 
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in demselben Intervall auch zwei Wurzeln hat. Wenn wir diesen 
Schluß mehrmals anwenden^ so finden wir, daß 93(0?) zwischen 
und 1 zwei Wurzeln hat. Das ist aber unmöglich, weil die Funktion 
3teii Grades (p^ix) nicht für vier verschiedene Werte von x (0 und 1 
inbegriffen) verschwinden kann. 

Es ist also bewiesen, daß die Gleichung ?P2*+i(^)^^ zwischen 

und 1 nur eine Wurzel x == — hat. 

Wir wollen jetzt die Gleichung 9?2a(^) "^ ^ untersuchen. Wenn 
diese Gleichung zwischen und 1 die Wurzel a hat, so folgt nach 
dem Rolleschen Satze, daß die Gleichung (P2ki^)^^ mindestens 
zwei Wurzeln zwischen und 1 hat. Das ist aber unmöglich, weil 

9^'2*(^) ^ 9^2*— 1(^) ^®* ^^^^ ^^ ®^®^ bewiesen haben, daß ^'g*— 1(^) 
nicht für zwei Werte zwischen und 1 verschwinden kann. 

Es besteht also der Satz: Für alle Werte von x zwischen 
und 1 hat die Funktion (p^j^ipo) dasselbe Vorzeichen. 

Um dieses Vorzeichen zu bestimmen, betrachten wir die Gleichung 

9^2*+i(^) ^ ^^ ^^® befriedigt wird für x==0 und x^-^' Nach dem 
Rolleschen Satze hat die Gleichung 

9^2*+i(^) = ^'ajfcC^) + ^2* = 
mindestens eine Wurzel a? = a, die zwischen und — liegt. 
Aus der Relation 

folgt, daß Ä2h nicht gleich Null sein kann, weil (f^^i^) nicht im 
Intervall (0, 1) verschwinden kann, und daß in demselben Intervall 
9?2^(a?) das entgegengesetzte Vorzeichen hat, als -4.2*. 

Es ist aber 

9?2* W — 1 . 2 . . . 2Ä; "^ -^1 1.2...(2Ä;-1) "*" ^ 1.2...(2Ä-2) 

Nach einem Satze der Algebra hat 9,-,(a;) für kleine positire 
Werte von x dasselbe Vorzeichen wie das Glied Atk—fr-n oder 

wie ^3jfe_2. 

Die Koeffizienten -4.2*— 2 imd J.2* haben also verschiedene Vor- 

I zeichen. Da aber A^== — >0 ist, so ist 



Zweites Kapitel. Die Jacob BemouUisclie Funktion. 45 

und für < a; < 1 besteht die Ungleichheit 

(-1)* ,,,,(«:) >0. 

Die Funktion 92*+i(^) ^^cl^solt ^^^ Vorzeichen im 
Intervall (0, 1), weil x = -^ ^^^ einfache Wurzel der Gleichung 

9^2*4-i(^) ^^ ^ ^®** ^^ diese Behauptung zu beweisen, nehmen wir 
an, daß 

sei. Die Funktion (p^^i^^) + -^2* verschwindet für x = a und x=^-^* 
Nach dem Roll eschen Satze soll 

verschwinden für einen Wert von x zwischen a und — ; was un- 
möglich ist. 

34- Bemoulllsche Zahlen. Die Koeffizienten A^,Ä^,A^, ,,, 
werden oft durch die Zahlen B^yB^yB^,,,. ausgedrückt, die man 
BernouUische Zahlen nennt. Es sind Zahlen, die mit dem Vor- 
zeichen + oder — , als Koeffizient von x in der ersten Potenz, in 
den Ausdrücken für 



>'^2 >^^ ^X% . . . 



auftreten. ^ o o 

Es ist 



^x'=1.2q>,(x) = l-2A^x^x^+^, 



X 

^ä:* = 1 . 2 . 3 . 4 g?6 (a?) = 1 • 2 • 3 . 4 ^4 a; + . . . 



X 

^ic® = 1 . 2 . 3 . 4 . 5 . 6 9?7 (ic) = 1 . 2 • 3 . 4 . 5 . 6 J^ a: + • • 



Es wird daher gesetzt 

1 2A^ = Bi, 1.2.3.4^4 = --B2, 1 2.34 5 6J« = -»3, . . . 

Allgemein ist 
oder . 1.2...2ä.^„ = (-1)-:B. 



(10) A,, = 



1-2-8...2«; 



Die ersten Bemoullischen Zahlen sind 

-«1=6' A=8Ö' -^»^42' -^*=P' -^8 = 66' 
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Die Tabelle der 62 ersten B findet man bei J. G. Adams im 
Journal für Mathematik, Bd. 85 (1878), S. 269—272. Die ersten 
20 Zahlen sind im Lehrbuch yon A.A. Markoff abgedruckt (S. 125). 

Auf die arithmetischen Eigenschaften der BernouUischen Zahlen 
gehen wir hier nicht ein * 

36. Entwlcklnng von 9>f^/^ix) In eine trigonometrlsohe 

Seihe. Aus der Theorie der trigonometrischen Brcihe (L. Dirichlet, 
Werke, Bd. I, S. 119 — 132) folgt, daß für 

0^a:<l 
die Entwicklung 

q>^j^(x) = "ö ö^o + «1 cos 2jtx + Og cos ^lscx -\ \- an cos 2nax -\ 

+ \ sin 2ücx + 62 sin 4xx -\ \-bn sin 2nücx H 

besteht. Die Koeffizienten haben die Werte 

1 1 

a„ = 2 / g?2ifc(^) ^^^ 2n%X'dXy bn='2 1 (p^^{x) sin 2n7cx • dx. 



Wenn wir x=l — y im Ausdruck för 6„ einsetzen, so wird 
wegen der Relation (8) in Nr. 32 



In 2 / g?j;t(y) sin 2n%y • dy. 





Pol^ch ist 6« == — 6n oder 6« = 0. 
Der Wert für 



»0 = ^ y8*(^)^^ 



läßt sich aus. der Relation . . 

bestimmen. Man findet 



X 

\ «0 = / VstC«)«^* Alk. 



* • 



... 

Um an zu berechnen, werden wir das Integral 

1 1 

/ 9^2* (^) ^^^ 2n%x ' dx durch / V^k-^^^ ^^® 2n^x • dx 


ausdrücken. Zu diesem Zweck gebrauchen wir die Formel 

* Diesem Gegenstand ist der dritte Abschnitt des Werkes von Herrn 
L. Saalschütz gewidmet (Yorlesnngen über die Bernonllischen Zahlen . . . 
Berlin 1893). 
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(11) jfftlf^'dx = g?^' — q)'t + ftlf(p"dx, 

die eine Folge der zweimaligen partiellen Integration ist. 
Setzt man in (11) 

(p = q)^j, (x), ^" = cos 2n7tXy 

so findet man für Ä > 1 

• / v ' , , sin2n9rri!; 



9>" == 9'2t-8(^) + ^2*-8, ^ = - 



2n9i; 
co8 2n9ra; 



Da aber 

1 



/ 



cos 2n7tx * dx^O 







ist, so wird 

1 1 





Wenn man diese Formel (h — 1) mal anwendet, so erhält man 
1 1 

/ 9^2* (^) ^^^ Swir'ic • dx = Bit—a / 9^2 (^) ^^^ Zwäoj • da;. 



Es sei jetzt in (11) 

(p = (p^{x)y ^" = cos 2n^x. 

Daraus folgt 

• 1 , • sin 2nyex 



9>"=1, ^ 



2 ^ 2WÄ 

cos 2n^x 



(2n«)> 

und 

1 



y 9?2(^) cos 2nücx • da; = ^g^- 





Man findet also 

(-l)*-^2 



^« "^ - n2* 



(2n5ry 

Wenn man iu der Entwicklung 

(p^Jix) *= "s- «0 + ^1 c^s 23ra; + öfg cos Atcx -\ 

a; gleich Null setzt, so findet man die bemerkenswerte Formel yon Euler 
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Diese Formel wird später benutzt, um die Konvergenz einiger 
Eeihen zu untersuchen. 

Aus den Werten für A^ und Ä^ folgen die Formeln 

■^ ^ 2« ^ 3* ^ 4* ^ 6 ' 

^ 2* ^ 3* ^ 4* ^ 90 

Dieses Kapitel habe ich unter dem Einfluß vonW.G. Imschenetzky 
(Kasaner Universitäts- Nachrichten 1870 und französisch in J. Hoüel, 
Cours de calcul infinitesimal, Paris 1888, Bd. 1, S. 476) und Herrn 
N. J. Sonin (Warschauer universitäts -Nachrichten 1888 und französisch 
im Journal für Mathematik, Bd. 116 (1896), S. 133) bearbeitet. 
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36. Bulersche Formel fOr ganze Funktionen. Es sei 

Wenn man diese Funktion yon bis x summiert, so findet man 



^f(x) =i)o9>„(^) + PlVn-l(p^) + • • • + Pn^i^ii^) + Pn-1^' 



Hier treten Bernoullische Funktionen ein. Wir wollen dieselben 
YoUstandig ausdrücken und die Glieder nach den Koeffizienten 
Ä^, A^, J-g, . . . anordnen. 

In dieser Weise entsteht die Formel 

^ ^ 



—8 ^n— 8 
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oder 



(1) 






Das ist die Eulersche SummationsformeL Sie ist genau fOr 
ganze Funktionen yom (n — 1)*^ oder niedrigeren Grade. 

37. Bestglied der BnlerMhen Formel. Wir wollen die 
Funktion 



(2) 



SP 



aosdrücken in dem Falle, daß f(x) keine ganze Funktion vom 
(n - 1)*«" Grade ist. 

Die Differenz von ^(a;) ist 



(3) 



X 

Afix) = f(x) - A lf{x)dx - AAf(x) - A,Af(x) 



^,_, A/"(— »)(«) - Än-iAf(^-'y(x). 

Aus der Int^ndrechnmig ist bekannt, daß 

AF{x) = Fix + l)-F(x)^F'(x) + :^F"{x) + :^^F"'(x) + ... 

1 





ist. 

Nachdem wir sukzessive 



X 

F(x) -fmdx, f^x\ fix), . . . , f(--^)(x) 



m=^n, n — 1, w — 2, . . . , 

gesetzt haben^ so yerwandelt sich (3) in 

Seliwanoff, Bififerenzenreclmaxig. 
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1 





1 



Am-^,m- ■ ■-T^^f^/'''-'K-)-ft^^/'K-^-)äu 





1 



-^f'\^)- ■ -r^:^/^-'K-)-ß^^E^/''K-+u)du 





1 



- A-if ("-^> (a;) -JAn-i{l - u)f<f^)(x+u)du 



oder 

Ai,ix) = - q>,{l)nx) - q>,(l)f"{x) -■■■ 

1 

- g',(i)f(-"(a') -Jv,a - «)fW(« + «)<^«- 



Da die Werte qPa(l), g?3(l), ..., 9>„(1) gleich Null sind, so ist 

1 



Nach der Formel (2) ist ^(0) = und folglich 

(4) V (x) = - / n (1 - «) 2/'^"^(* + «)<^«*- 



Es wird hier nach x summiert nnd nach u integriert 
Wir werden n = 2Ä setzen* Wegen der Relation 

ist dann 

(5) ^(«) - - ffp^^(u)du'^f('''\x + «). 

Das ist das Bestglied Ton Jacobi (Journal fOr MatL, Bd. 12 (1834), 
p. 263 oder Werke, Bd. 6, p. 64). 



Drittes Kapitel, Enlersche Stunmationsformel. 5I 

Wenn man in (5) die Bernonllisclie Funktion (p^jj^^ durch 
die trigonometrisclie Reihe (Nr. 35) ausdrückt^ so findet man das 
Restglied von Poisson (Pariser Memoires, Bd. 6 (1826), p. 580). 
Wir werden die Jaco bische Form gebrauchen. 

Die Funktion 9?,^(m) behalt ihr Vorzeichen im Intervall zwischen 
und 1. Deswegen folgt aus (5) 

1 



wo < e < 1 ist. 
Da aber 

so ist 



und folglich 



1 




Je 

(6) ^(a;) = ^,*2/^")(a; + 0). 



Wir haben also bewiesen^ daß ftir jede Funktion f{pc)y die 
gewissen Stetigkeitsbedingungen genügt, die Formel besteht 



(7) 



m + /•(!) + /•(2) + . . . + f{x - 1) 



« 



(8) 



=Jf(x)dx + Wf{x) - f(p)] +A,[f'ix) - f'(Oy\ + . • ■ 

1+ ^2*[/'("'(e) + /•(")(e + 1) + . . . + f(»>')(B + x- 1)]. 

Für a; — 1 i8t 
1 

m ^ff{x)dx + Mm - rm + Mf'(}) - fm + • • • 



A,i-i[p*-»)(l) - /•(»*-»)(0)] 



Die Qlieder mit Än—i fallen weg, sobald % > 1 ist (Nr. 32, 
Form. 9), Für * = 1 ist 

4* 
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m+f(i)+f(2)+..+f(x-i) 



(9) • 



X 



^ff{x)dx+A\f{x)-m-\ 



und 



\+\\f'iP) + f'(ß + 1) + f"(ß + 2) + . . . + r(ö + « - 1)] 



(10) f(0) =J*f{x)dx + A,[f(l) - fm + Af'(ß)' 



88. Andre Form des Bestgliedes. Wenn die Ableitungen 
/•(")(<) und /•(«*+ «)(0 für Werte von t zwischen und a; ihre Vor- 
zeichen behalten und dieselben Vorzeichen haben^ kann das Rest- 
glied (6) eine andre Form annehmen. 

In der Formel 

/ q>o"du == (po' — y'cD + / oq>^'du 
setzen wir _ 

9 ^^f^'^Kx + m), a" = 9>„(m) + Au. 





Dann wird 

X 



y'=2/^"+«(a; + «), o' =?.„+,(«), 





X 



und folglich ist 



-f^P'K«' + «) • [»»„W + ^tkldu 



« 








. Es besteht also die Relation 

1 





1 



-jvi.i«)^^ •2/'^"^(« + «) 
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Da aber 
1 





80 ist 



/' 



fm(x + u)du - /•("-«(« + 1) - /^«*-i)(a!) =. A /•("-«(«), 



tß 



und folglich 

1 



V 



) 



Das erste und letzte Glied dieser Formel haben verschiedene 
Vorzeichen, weil, nach dem in Nr. 33 bewiesenen Satze, ^>^^{y) nnd 
9^2*+sW ™ IntervaUe (0, 1) verschiedene Vorzeichen haben und 
nach unsrer Voraussetzung f^^*^\t) und /*^**~^)(^) gleiche Vorzeichen 
haben. 

Die gefundene Relation hat die Form 

a^b + c, 
wo ac < ist. Aus der Qleichung 

a^^ ab + ac 

folgt 

0<a«<a6, 0<J<1. 

Es ist also A, 

a=^db, 

wo < e < 1 ist. 

Nach dieser Bemerkung ist 

(11) ^(o?) = e . -42*[/^"-^K^) - /•^"-^>(o)]- 

Diese Formel kann statt (6) gebraucht werden, sobald f(t) gewissen 
Bedingungen genügt (Vgl. Markoff, Differenzenrechnung, S. 121.) 

89. Zweite Form der Bulereohen Formel. Wir werden 
voraussetzen, daß f&r alle positiven Werte von t die Ableitungen 
f\t), /'^^(O, f^^{t), . . . dieselben Vorzeichen behalten und daß aUe 
Ableitungen f(t), f"(t), f"\t), . , . unendlich klein werden, sobald t 
unendlich groß wird. 
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In diesem Falle hat das Integral 

1 

00 








einen Sinn. Es bedeutet, daß das Integral 

1 








bei unendlich wachsendem x sich einer bestimmten Grenze nähert. 

In der Tat, bei wachsendem x vergrößert sich der absolute 
Betrag des Integrals, weil neue Elemente mit demselben Vorzeichen 
zutreten. Sobald aber x hinreichend groß ist, so ist 

1 



•/ 



<|^2*|{l+|/'("-^>(0)|}, 



weil hier die Formel (11) anwendbar ist und nach der gemachten 

Voraussetzung |/"(**-^)(a;) | beUebig klein sein kann. 

Nach bekanntem Prinzip hat das Integral eine bestimmte 

Gtrenze. 

Aus der Relation 

1 

1 
folgt « 



«/ 





1 



•^ .00 



Da aber 
1 



/» W y^ OO 

<p,,(.u)du^fi^^)(t + u) und / v,,+Md»y!P'+'Ki + «) 

Terschiedene Vorzeichen haben, so ist nach der Bemerkung im 
vorigen Paragraphen 



• ' ,•» * 
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r<p,,(u)du^f(»'')(t + «) = 01 ^«/■("-i)(o), 

« 

wo < 01 < 1 ist. 

um diese Formel zu Terallgemeinem, genügt es zn setzen 

m = /(^ + «) = j(«), 

wo t; eine neue unabhängige Variable ist und x einen bestimmten 

Wert hat. 

Nach der bewiesenen Formel ist 
1 

(p,,{u)du^F^^^v + m) = 02 ^2*-F("-^K0) 



oder 

1 

* 

wo < 02 < 1 ist. 

Wenn man in der Eulerschen Formel das Bestglied 

1 





durch 







OO j<^ CO 



•^ X 





ersetzt; so findet man 

' ~ + 



X 

•/ 



» 



• • • + ^„_i /'("-»)(«) + / iPiti»)d^2f^"Kt + «)» 



wo 

1 

" 

Jetzt woUen wir beweisen; daß die Eonstante Ck von Je un- 
abhängig ist. 
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Für * = 1 ist 

X 



2/"(0 - fmdt + C, + AJ(x) + ß,(u)du^f"(t + «). 





Mit Benutzung von (12) findet man 

y!m = c» + imdt + Ä^fix) + A^f\x) + • • • 

tind , • 

» sr 

Der Vergleich der beiden Formeln gibt die Relation 

Ci + A^f (x) -{■ ■ ■ + Ätu-if(*^-%x) + e ^,i /■("-»)(«) 

^C, + B^A^f{x), 

die für jeden ganzzahligen positiven Wert von x besteht. 

Bei nnendUch wachsendem x sind aUe Ableitungen von f{x) 
unendlich klein und folglich ist Ck ^ C^. 

Es ist also die Formel bewiesen 

X 

^m=c+fmdt+Ä,f(x)+A^f'(x)+- 

wo < ö < 1 ist. 

40. Der allgemeiner« Fall der Bnlersohen FormeL Bis 

jetzt haben wir die bestimmte Summe von der Form 

f(0) + f(l) + f(2)+... + fix-l) 

betrachtet. Nun wollen wir zur Summe 

F(a) + F(a + h) + F(a + 2h) +-"+ F(b - h), 

wo b == a + xh ist, übergehen. 

E« ^' F{a + th) - m; 

dann ist 

FQ> - A) = f{x - 1) 
nnd folglich 

F(a) + F{a + h) +...^ F(J>-h)=f{Q)) + fO.)+-- + fix-l). 



(13) 
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Wenn wir beracksichtigen, daß 

/•(«)(<) = Ä" J'W(o + th), 

80 nehmen die Formeln (7), (8), (9) nnd (10) die Formen an 
■F(a) + F(a + h) + F(a + 2h)-\ i-Fib — K) 
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(14) 



- 



^/^F(u)du+Äi[Fih)-Fia)]+J^h[F'(h)-F'(ay\ + - ■ ■ 



(15) 



+ .l,t_i Ä»»-» [J'("-»)(6) - 2!'(»*-«)(a)] 

+ ^,tÄ»[J'(»)(o+eA)+J'(»*)(o+l+eA)+- ■■+F**(h-h+eh)l 

a+k 

F(a)^^ I^F{u)du+J^[F(a+h)-F(a)]+J^h[F'{a+h)-F'(ay]+- • • 

+i,i_iÄ»*-»[2!'<»*-»)(o+A)-^<»»-«)(a)]+ui«Ä«»Z{»»)(o+eÄ), 
.F(o) + Z(a + Ä) + Z(o + 2A) + • • • + J?'(& - Ä) 



(16) 



jfF(u)du +Ä, [F(b) - Fia)] 



und 



+ J,Ä» [J"'(a+eÄ) + J"'(a+1 + 0Ä) + •• • +F"(b-h+Bh)] 

a+h 

(17) 2!'(o) = ^ JF(u)du + ^1 [2!'(a + A) - i?'(a)] + J, A*2!"'(a + 0Ä). 

a 

In alleil diesen Formeln ist & » a + a;A' nnd x eine ganze 
positive Zahl. 

Wenn die Ableitungen F"{{), F^(t), JPW(<), ... für ^^a die- 
selben Vorzeichen behalten und die Ableitungen F*(t\ F"(t), F'^^t),. . . 
unendlich klein werden, sobald t unendlich groß wird, daÄn besteht 
die Formel 



(18) 



F(a) + F(a + h) + F(a + 2h) + - ■ ■ + F(b-h) 

b 

'=C+^rF(u)du+AF{b)+J^hF'(p) + --' 
+ Ä,t-i Ä"-» J'(»*-»)(6) + ö^,t A»*-»l?'("-y(6), 



wo die Konstante G yon b nnd h unabhängig ist. (Vgl. Markoff, 
Differenzenrechnnng, S. 131 — 133.) 
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Viertes Kapitel. 
Anwendungen der Eulerschen Formel. 

41. Entwicklung von — ^ — In eine Belhe. Setzt man in 

e^— 1 

der Formel (8) (Nr. 37) 

so wird 

1 



/ 



e*-t 



ßt'dX'^-^, /■«(!) = f'»«', /'W(0) = ^'», /•(")(«!) = ^«*c«« 





und folglich 



e*-l 



^ +A(c'- l)+^<(e'- !) + ••• 
■ ■ Nach der Formel (12) (Nr. 35) ist das Bestglied 

Bei wachsendem Je bleibt 

endlich, weil 

oder 1 

i<M,< — ^— 

1- 



22*-l 



ist. Daraus folgt, daß lim ilf^ = 1 ist. 

Es ist also für |^|<27i; das Restglied dem absoluten Betrage 
nach beliebig klein für hinreichend große Werte von Je. 

Folglich besteht für |^| <23i; die konvergente Reihe 



e'-l 



1 = ^ + A (e' - 1) + AKe' - 1) + At\e* - 1) + • • • 
oder 

(1) -^ = }+^, + ^^ + ^,^» + ... 

42. Entwicklung von ctga; und tga; In eine Belhe. 

In der Formel {8) (Nr. 37) setzen wir 

Dann wird ^^^^ = '''' (^^> 



1 



/ 



cos (tx)dx = ^^, f (t) t sin (tx) , f" (0 = + <* sin (tx) , 

. . ., /•(»*)(a;) = (- If-H^" cos {tx), 
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und folglich ist 

1 =• ?^ + Ä^ (cos < — 1) — A^t sin ^ + AJ^ sin ^ 

+ Äik^2(- l)*-i<2*-8 sin* + ^2,(- ly-H^^ cos (ßt). 
Da aber 

1 - J.i(cos ^ - 1) = -^ cos« -, 

so findet man nach der Division durch 

sin ^ = 2 sin 2^ cos -^ 
die Formel 

cos (9 *) 



+ (- l)*-i^,*_,<"-»+ (- 1)*-M,i<" 



sin^ 



Ebenso y wie in Nr. 41, kann man beweisen, daß das Bestglied 
für 1^1 < 2ä und für hinreichend große Je beliebig klein wird. 

Wenn man x statt -^ einsetzt, so entsteht die Entwicklung 

(2) ctg o; = - - 2^A^x + 2^A^x^ - 2^A^x^ + • • •> 

sc 

die für '\x\ <.% konvergiert. 
Mit Hilfe der Relation 

tg a; = ctg a? — 2 ctg 2x 
findet man die Reihe 

(3) tga; = 2«(2^-l)^ir;-2^(2*-l)^^a?8+2«(26~l)^a;'^ ; 

wenn | a? 1 < ^ ist. 

48. StLrlingsche Beihe. Wir wenden die Formel (18) (Nr. 40) 
an und setzen 

l^(ri;) = log rc, a=l, ä = 1, b = x. 
Man findet dann 

/ log xdx = xlogx — x, F\x) = x-^y F^\x) = — xr-^, • • • 

..., 2^H(^)==(~l)"»-U.2.3...(m-l)Är-"» 
und folglich 

log [1 • 2 • 3 . . . (rc — 1)] = C + a; log a; — ic + J-i log x 

+ J^ar-i + ^^1.2ar-8 + --- + ^2*-srl-2-3...(2Ä-4)Är-(«*-«) 

+ 6^2*1-2. 3. ..(2Z;-2)a;-(2*--i). 
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Indem man zu beiden Seiten log x addiert nnd Ä^ durch 
seinen Wert ( — -r-j ersetzt, so erMlt man 

log (1-2-3 ...«) = C+ (« + i) log a; - a; + J,«-* 

(4) I + ^J • 2«-» + • . • + ^,t_,l • 2 • 3 . . . (2* - 4)a;-(«*-») 

+ e^jt 1 • 2 • 3 . . . (2* - 2)a;-(»»-»). 

Hier ist J; > 1 vorausgesetzt. Für k = l besteht die Formel 

(5) log (1 .2-3 ... «) = C+ (a; + I) log a; -T- a; + %^A^x-K 

• • • 

Es sind und 0^ Qrößen zwischen und 1. 
Es bleibt noch die Eonstante C zu bestimmen. Zu diesem 
Zweck wenden wir die Wallissche Formel an: 



2 



= lim ^(n) für lim » = oo, 



wo 



-/ V 2-2 4-4 6-6 2n-2n 

.. ^ >' l-S 3.6 6.7 (2n-l)(2n + l) 



ist Mittels ein&cher Umformungen erhalten wir 

ih( \_ r 2.4.6...2n " i« 1 
W "~ Li • 3 . 6 . . . (2n-l)J * 2n + 1 

_ r2 ■ 2 . 4 • 4 • 6 • 6 . . . 2n • "^n ^ 1 
""L l-2.3.4...(2n-l)2n J *2n + l 

(2.4.6...2n)* 1 2*". (1.2. 8... n)* 1 



(1.2.8.4...2n)* 2n + l (1 • 2 • 3. . .2n)*.2n , , 1 

Daraus folgt 
log *(w) == 4n log 2 + 4 log (1 • 2 • 3 . . . ») 

-21og(1.2.3...2»)-log(2n)-log(l + ^) 

oder mit Benutzung von (5) 

log*(w) = 4n log 2 + 40 + (4» + 2) log» - 4n + ^^lA^^ 

- log 2 - logn -log(l + ^) 

- (4» + 1) log 2 - 2(7- (4n + l)'logn + 4» - e,^,-^ 
= -21og2 + 2C+4e,J,i-log(l + ^)-0,J,i. 
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Wenn wir jetzt n unendlich wachsen lassen, so wird 
lim log *.(n) = log |; lim ^4d^ J, -J =: 0, 

nnd folglich 

log| = -21og2 + 2a 

Daraus findet man 

C=«logy2Ä. 

Es bestehen also die Formehi 



(6) 



log (1 . 2 . 3 ... a;) = log y^ + (x + ^)log x - x + A^x-^ 

+ Ä^l'2x-'^ + " +-4,*«,1.2.3...(2Ä-4>-(«*-»> 

+ 6-42*1 •2-3... (2* -2>-("-^) 
nnd 

(7) log(1.2.3...a;) = logi/23r + (a; + |)loga?-a;+6i^a;-^ 



1 . 



Da < 01 < 1 nnd -4s, = ^ ist, so folgt ans (7) 

log 1/2^ + (a; + -^) log a; — a? < log (1 • 2 . 3 . . . a;) < log y2x 

oder • ^ ^/ . ,..x-^ . . ^ 

(8) y2jtx • a^ • 6-* < 1 • 2 • 3 . . . a; < Y^tcx • ir* • 6-* • 6"*. 
Für große Werte von x kann die Nahernngsformel 

(9) l-2-3...rp=->/25r^.rr*.e-* 
gebraucht werden. 

44. Untersuchiing des Sestgliedes der StLrlingschen 
Beihe. Die Gleichung (6) ist die Stirlingsche Reihe, deren Eon- 
yergenz wir jetzt untersuchen 'wollen. Zu diesem Zwecke muß man 
sehen, ob das Bestglied bei unendlich wachsendem Je unendlich 
klein wird. 

Nach der Formel (12) (Nr. 35) ist 

BA^jt •1-2 . 3 .V. (2* - 2)rp-("-i) 
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In Nr. 41 ist bewiesen, daß bei wachsendem Je sieb die Reibe 

der Grenze Eins nabert. Für jeden Wert von x wird aber das Produkt 

12 3 2Ä;-2 



• • • 



27CX 2nx 27CX 2%x 

beliebig groß für binreicbend große Werte von Tc, 

Daber ist die Stirlingscbe Reibe divergent. Für große Werte 
von X und kleine Werte von Je ist das Restglied dem absoluten 
Betrage naeb sebr klein. In diesem Falle gibt die Stirlingscbe 
Reibe eine gute Annaberung für 

log(l-2.3...Ä;> 

Z. B. für a; = 100 und Z; = 3 ist der Febler < — ^- (Markoff, 
Differenzenrecbnung, S. 137.) 

Für einen bestimmten Wert von x ist der Febler dem absoluten 
Betrage nacb kleiner als 

|^2*|-l-2.3...(2Ä;-2)a;-("-i). 

Die beste Annaberung bekommt man für einen solchen Wert 
von Jcj für den der Ausdruck 

— 9 MV 2 8 2fe-2 fi , J_ , J_ , 1 

den kleinsten Wert annimmt^ weil der absolute Betrag des Restgliedes 
kleiner als Uk isi 
Für Ä ^ 6 ist 



22* • 32* ' ^2* • 58* • ^ 2047 

Deswegen kann man angenähert setzen 

o / 1 \8 2 3 2Ä;-2 

X^icxj 2'xx 2nx 2nx 

Da aber angenähert 

1 \« 2 2Ä;-2 2h-l 2k 



^*+^='2a;.(^).^ 



• • • 



2nx 2nx 2nx 
SO ist 



^ifc+i ^ 2fe-l ^ _2k_ /k\^ 



u^ 27CX 2 



«05 \7CX/ 
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Die absoluten Beträge der Glieder der Stirlingschen Reihe 
nehmen also ab^ solange h Kjcx bleibt. Die beste Annäherung erfolgt 
also, wenn k nahe zxi %x + 1 liegt.* 



45. Angen&herte Berechnung der Summe 

1 + 2 + 8H ^T=^ 

ffir groCbe Werte von x. In der Formel (18) (Nr. 40) setzen wir 

F{x) — —} a=lj Ä = 1, h = x. 

Dann finden wir 
dx 



ß 



= log X, F^'^\x) = (- 1)«1 . 2 • 3 . . . Wir-<"»+i) 

und folglich 

1 + 4 + 1 + ••• + ^ = C + logx-ii-Aa.-» 

-A^12- 3a;-* ^2i-jl • 2 • 3 . . . (2ä - 3)a;-("-») 

-e^2il-2-3...(2Ä-l)a;-»* 



(10) 



Die Eonstante C läßt sich nicht durch bekannte Eonstanten aus- 
drücken, sie läßt sich nur angenähert berechnen. 

Wenn man o; =» 10, Ä; = 7 setzt und das Restglied vernachlässigt, 
so findet man 

C- 0,57 721 56649 01532... 

Das ist die Eulersche Eonstante. 

Ebenso wie in Nr. 43 kann man zeigen, daß die Reihe (10) für 
jeden Wert von x divergiert. Deswegen darf zur Berechnung der 
betrachteten Summe nur eine geringe Anzahl der Glieder der Reihe 
genommen werden. 

In die Untersuchung, wie groß diese Anzahl sein kann, werden 
wir uns nicht einlassen. 



• Vgl. Encykl. der Math. WisB., Bd.I, p. 103. 
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Differenzengleiclinngeii. 



Erstes Kapitel. 

Allgemeines filier Differenzengleichimgen. 

46. Zwei Formen der Dlfferenzengleloliimgen. Eine 
Oleichnng von der Form 

(1) *(^; y^y ^y^y ^^y^y • . • ; ^"•yj « 

heißt Differenzengleichang. Hier ist $ eine gegebene Funktion, 
y eine unbekannte Funktion. 
Bei Anwendung der Formeln 

^*y, - Vx+t - 2y.+i + y,» 



nimmt die oleichnng (1) die andre Form an 

(2) ti>(x, y^, y,^„ y,^„ . . . , y,^„) = 0. 

Ist die Gleichung in der Form (2) gegeben, so kann man sie in 
die Form (1) überfahren. Es genügt, die Formeln 

y«+« = ». + »» Ay, + *"^^ ~ ^^ A»y, + • • • + A-y, 

anzuwenden. 

Von diesen beiden Formen (1) und (2) ist (2) die bequemere. 
Es ist einfacher, eine gegebene Funktion statt y in die Gleichung (2), 
als in (1) einzusetzen. 
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Wenn bei der Transformation (1) in (2) z. B. die Größen y^ 
und j/^ , ^ wegfallen^ so wird die Gleichung (2) dadurch vereinfacht, 
daß man y , «==;5f setzt. 

Es sei die Gleidmng 

^\ + A\ - Ay^ - y^ = 
gegeben. Mit Anwendong der Formeln 

geht diese Gleichnng in 

über. Setzt man y, ■ , = ^x' ^^ entsteht die einfachere Gleichung 

Aus diesen Gründen werden wir voraussetzen, daß die gegebene 
Gleichung die Form (2) hat und daß in ihr y^ enthalten ist. Wenn 
außer y^ noch ^^^ i^, aber kein y mit größerem Index vorkommt, so 
werden wir (2) eine Gleichung m*®' Ordnung nennen. Die eben 
betrachtete Gleichung 

A%+A\--Ay^-y^:^0 
ist also eine Gleichung 1*®' Ordnung. 

47. Parttkul&re und allgemeine Lösungen. Wir werden 
voraussetzen, daß die vorgelegte Differenzengleichung w*®' Ordnung 
die Form hat 

Die Funktion f soll für alle betrachteten Werte von x, y^, Vx-i-v 
3/ar+2^**' yx+n—i ®^®^ endlichen und bestimmten Wert erhalten. 

Jeder Ausdruck y^^ ^{pc)^ der dieser Gleichung genügt, heißt 
partikuläre Lösung von (3). 

Es gibt unendlich viele partikuläre Lösungen, um dies zu zeigen, 
wollen wir den Größen j/q, j/^, y^y-j yn—i '^ilU^^li^^liö Werte 
beilegen. 

Aus (3) ergibt sich für a; = 0, 1, 2, 3, . . . das System der 
Gleichungen 

Seliwanoff, Differenzenreohnung. 5 
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W iyn+8=/*(2^ y^y y^y y^y'"7 Vn+O^ 

yn+z^fi^y Vsy y^y %; •-; »«+2); 

aus welchem y^, y„ 1 1; ^»+2^ ^«+8^ • • • sukzessive berechnet werden 
können. Es hat aliso y^ einen ganz bestimmten Wert für ganze 
positive X, sobald j^^; y^; y^y ---y ^n— 1 g^g^^en sind. Man kazm also 
schreiben 

(5) y,= F(Xy yoy yiy y^y-'y y„-i)- 

Die Wahl der Anfangswerte y^, y^, y2;«**; y«— 1 ist nur einer 
Bedingung unterworfen: es sollen alle Gleichungen des Systems (4) 
einen Sinn haben. 

Es ist möglich; daß die Funktion 

f{^y yxy y«+i; y^^iy-y y»+«-i) 

einen bestimmten Wert erhält fQr rc^a und für beliebige Werte 
von y^, y^^i, y^^,, ..., y«+„.r 
Dami besteht das System 



(6) • 






Daraus folgt die Relation 

fOr a; ^ a. Hier ist die Differenz x — a gleich einer ganzen posi- 
tiven Zahl. 

In den Formeln (5) und (7) sind y^, y^, y^, ..., y„_i oder 

y«' Va+v ya+v ••■> y«+«-i willkarUche Konstanten. 
Jede Lösung der Gleichung (3) Ton der Form 

(8) y.= ^(«, c?i, q,, C„ .... C,), 

die n willkürliche Eonstanten enthalt, heißt allgemeine Lösung, 
wenn es möglich ist, die Gleichung (8) in die Form (5) oder (7) 
zu bringen. 
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Zn diesem Zweck mnß man aus dem System 



y„_i=^(*»-l; Ci, Ca. • • v C;,> . 



(9) oder (10) 









^y«+n-i=^(«+W-l; C,, C^, . . ., C„) 



die Konstanten C^, C^, . . ., C^ durch y^, ^i? ^2; • • •; y„_i ^^^^ durch 
ya> ^a+i^ y«4-2^ • • -^ ^a+n-i »^sdrücken. 

Man nennt die willkürlichen Konstanten C^, C^, - - -, Cn wesent- 
lich, wenn das System (9) oder (10) in bezug auf diese Konstanten 
lösbar ist. 

Allgemeine Lösung von (3) ist also eine Lösung^ die n 
wesentliche Konstanten enthält. 

Aus einer allgemeinen Lösung lassen sich viele andere ableiten. 
Es genügt in (8) 

(Ä= 1, 2, . . ., n) 
zu setzen. Die Koeffizienten a^^ sind nur einer Bedingung unter- 



worfen, daß die Determinante 



a 



kh 



(Ä; = 1, 2, . . . , n; Ä = 1, 2, . . . , w) 



nicht verschwindet. Li diesem Falle lassen sich c^, c^, . . . , c„ durch 
(^v (^29 ' - -f C7« ausdrücken. 

48. EigenBChaften der linearen Dlfferenzengleichnngen. 

Eine Gleichung von der Form 

(11) y.+n+Piy.+n-i+p^yx+n-2+"'+Pn^iyx+i+Pnyx= Q. 

heißt lineare Gleichung. Die Koeffizienten p^, 1>2>---? P — 1> P 
und das letzte Glied Q^ sind gegebene Funktionen von x. Wenn 
wir Q^ durch Null ersetzen, so entsteht die entsprechende 
homogene Gleichung 

(12) y.+« + i>iy.+n-i + Ay«+n-2 + • • • + Pn-iy.+i + Pny^ = 0. 

Die Gleichung (11) heißt eine vollständige, wenn Q^ nicht 
gleich Null ist. 

Es sei Yx eine partikulare Lösung von (11). Wenn man in 
die linke Seite dieser Gleichung 

(13) y.-y,+^. 

einsetzt, so findet man 

5* 
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Die erste Zeile ist nach der Yoraussetzung gleich Q^. Daraus 
folgt^ daß der Ausdruck (13) der yollständigen Gleichung genügt^ 
wenn 0x ©ine Losung der entsprechenden homogenen Gleichung (12) ist. 

allgemeine Lösung von (12) ist; dann ist 

yar= ^x+ V{^y ^V C%^ . • •; Gx) 

aUgemeine Lösung der vollständigen Gleichung (11). 

Lineare homogene Gleichungen besitzen folgende bemerkenswerte 
Eigenschaften: 

1. Wenn z^ eine Lösung von (12) ist, so ist auch C^^ 
eine Lösung. Dies folgt aus der Relation 

Hier ist G eine willkürliche Konstante. 

2. Wenn z^^ und z^'^ Lösungen von (12) sind, so ist auch 
;8f(i) -j- ;8f(2) eijie Lösung, weil 

+ Pn(4'' + 4'') = {4%n + Pl4%n-t + Pi4%n-, + -+ Pn'^f') 

ist. 

Aus diesen Sätzen folgt der wichtige Satz : 

3. Sind z^\ 8^\ ..., z^^ verschiedene Lösungen der 
linearen homogenen Gleichung, so ist auch 

(14) y^^ C,gm + C,^W + • • • + C 4«) 

eine Lösung. Hier sind C^ C^, . . ., C« willkürliche Konstanten. 
Wenn dieselben wesentlich verschieden sind, dann ist der Ausdruck 
(14) aUgemeine Lösung von (12). 

Li Nr. 47 haben wir gesehen, daß eine Differenzengleichung eine 
Funktion definiert, wenn die Anfangswerte gegeben sind. Um n 
partikuläre Lösungen zu bilden, genügt es, n Systeme von Anfangs- 
werten anzugeben 



S: 



(15) 
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^0 ' 1 > ^^2 > ' ' ' f '^n — lf 
Z^^\ ^W, i?(2) J^(2) , 
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;^(«), j^(«) ;^(W) ;^(«) . 



Wenn die Determinante^ die ans diesen Elementen gebildet ist^ 
nicht verschwindet; dann sind im Ausdruck (14) die Konstanten 
wesentlich verschieden und die Formel (14) ist wirklich allgemeine 
Lösung der Gleichung (12). 

Wir bekommen also den Satz: 

Die allgemeine Lösung einer linearen homogenen 
Gleichung ist eine lineare homogene Punktion der will- 
kürlichen Konstanten. 

Es war bei dem Beweise dieses Satzes vorausgesetzt; daß die 
Koeffizienten der gegebenen homogenen Gleichung einen Sinn haben 
für ganze positive x. 

Wenn aber wir nur behaupten können^ daß i>i; JP2; • • •> P bild- 
liche und bestimmte Werte haben für x'^a, dann ist das System (15) 
durch das folgende 

a ' a+1' ' a-\-n — 1' 



(16) 









ZU ersetzen. 

Nach diesen allgemeinen Betrachtungen gehen wir jetzt 
speziellen Gleichungen über. 



zu 
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49. Homogene Oleichung. Um die aUgemeine Lösung der 
Gleichung 

zu bilden^ geben wir der Größe y^ einen beliebigen Wert. Aus 
(1) folgt dann 



y«"" A-^l-^ • • • -^ar-l^O- 
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Weil ^0 willkürlich ist^ kann man schreiben V 
(2) y,- C{A,'Ä^'A, . . . A^^^^ ^' 

wo C die willkürliche Konstante ist. Das ist die gesuchte Lösung. 

Wir werden Toranssetzen^ daß in (2) keine der Ghrößen A^^ A^^ 
A^y ..., Ax^\ gleich Null isi Sonst hätten wir die eyidente 
Lösung y^= 0, die kein Literesse bietet. 

Wenn für alle x^a die Funktion Ax einen endlichen; be- 
stimmten und von Null verschiedenen Wert hat; dann ist die 
allgemeine Lösung unserer Gleichung 

(3) y^= C * Aa^ AaJ^X' AaJ^% . . . Ax^\> 

50. Vollständige Oleichnng. um die Gleichung 

(4) y.-^x-Ky.-^^ 

aufzulösen; setzen wir 

(5) y,,= Wx-v*, 

wo Ux und 'Ox neue unbekannte Funktionen sind; eine derselben 
wählen wir sO; daß die Gleichung für die andere Funktion möglich 
einfach sei. 
Da 

t?«+l == ^* + ^ ^'ar 

ist; SO nimmt die Gleichung 

UxJ^\ • Var+l — AxUx'Ox == ^x 

die Form an 

Waj+l(t?ar + Av^r) — Ax^x'Ox == Bx 

oder 

Wenn Ux eine partikuläre Lösung der entsprechenden homogenen 
Gleichung (1) ist; dann erhalten wir die Gleichung 

die Vx bestimmt. Es ist 
und folglich 

Aus der Betrachtung in Nr. 49 folgt; daß man 

Ux = A^A^A^ . . . Ax^x oder Ux = AaAaAf.\ . . . Ax~~\ 
setzen kann. Li diesen partikularen Lösungen ist Uq = 1 oder Ua ^ 1. 



\. 
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51. Beispiele. 1. Es sei die Gleichung 
gegeben. Hier ist Ax^5, JB« = a*, 



Die Formel 



^'-y.7Sr. + 



^8*+ 



hat verschiedene Werte, je nachdem a = 3 oder a ^ 3. 
Wenn a Ton 3 yerschieden ist, dann ist 

3 

und folglich 

Wenn aber a =-= 3 ist, so ist 

nnd t 

2. Wenn wir die Gleichnng 

, Sx-{-l X 

Vx^i "T- 3^-4-7 ^« ~ (8a; + 4) (8a; +7) 
nehmen]^ so ist 

^^ 8ic + 7' ^* (8a; + 4)(8a;+7/ 

Wx — (, ^A • Y • 10 " 13 ' 16 * " 3a; 4-1 ' 3a;4-4 



= (- 1> 



1+1 30! + 

1-4 



(8a! + l)(8a! + 4) 

i = _i(_l)..a;, 
««+1 *^ '^ ' 



27(- 1)« -^ = 2?a5 . A (^ = - ia;(- 1)-- i2;(- 1> 

Folglich ist die allgemeine Lösung dieser Gleichui^ 

^'^ ("" ^^ * (3a; +1) (3a; + 4) ' ^■'" \2^ " ij ' (8a; + 1) (8a; + 4) 



^ 

1 
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3. Als letztes Beispiel wollen wir noch die Gleichung 

20? + 7 Sg + l __ 1 

auflösen. Es ist 

7 9 11 2x + 6 
'^^ 1 3 5 2a;- 1 

14 7 3a;-2 (2a; + 1) (2a; + 3) (2a; + 5) 



• • • 



4 7 10 Sas+l 16(8«+!) 

^, 15 

«♦«+1 (2« + 8)(2a + 5)(2a; + 7)' 

'ei 1 i'^ L 

^(2a! + 3)(2a; + 5)(2a! + 7)'" 8^ 
1 1 



(.+|)KI) 



+ c^ — - - h c. 



10/ I ^\/ I ö\ * (2a; + 3)(2a; + 6) 



^ (2a;+l)(2a; + 3)(2a; + 6) p_£ 2a; + l 
^«"" 3a;4-l 4 * 3a;+l' 

(Vgl Markoff, A.A. Differenzenrechnung, S. 156 — 158.) 
Jetzt wollen wir die Anwendungen der Gleichungen erster 
Ordnung zeigen. Die Koeffizienten yerschiedener Entwicklungen 
genügen Differenzengleichungen. Die Auflösung derselben führt 
zur Bestimmung der Koeffizienten. Als Beispiel nehmen wir die 
binomische Entwicklung und die Entwicklung von cos xt nach Potenzen 
von cos t 

52. Die blnomisohe Bntwickliing. Das Produkt von x Faktoren 

hat die Form d + (1 + • • • (1 + 

Der letzte Koeffizient ist gleich Eins. Die Funktion 

Ax hat einen Sinn für a? ^ 1, 

•^x V yf t> f) ^ ^ ^} 

^X )9 }} ff yf *^ ^ ^9 



Die Anfangswerte sind 

-4i=l, £2=1, Cgsal,... 

Nach Multiplikation mit 1 + ^ findet man 
(1 + ^)«+i= 1 + (1 + A^jt + (A,+ B,)fi+(B,+ C,)fi + 
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Anderseits ist 

(1 + ty^i == 1 + A^^tt + B,+i^ + C^j^tf + . • • 

Daraus folgt 

Gc+l = Caj + Ba-, 



Die Auflösung dieser Gleichungen gibt 

Ax= 21 = a? + c; es ist c = 0, weil -4^ == 1 ist. 

Aus der Bedingung ^Bg = 1 folgt c = und 

Die willkürliche Konstante c ist = 0, weil Cg = 1 ist. 
In dieser Weise erhalten wir die binomische Entwicklung 

... (i + ty=.i + xt + ^J^fi+^J^^^ 

x{x-'l){x-2),..[x-{X'-i)] 
■*" 12 3...« 

Der letzte Koeffizient ist gleich Eins^ wie es auch sein sollte. 

58. Entwicklung von coa x t nach Potenzen von coa t. 

Nach der Mo ivr eschen Formel ist 

(cos < + * sin ^)* = cos xt + i sin xt 

Der Vergleich der reellen Teile beider Seiten dieser Gleichung 
ergibt 

cos xt = cos*^ — ^^f "I ^ sin^^ • cos*-"^^ 

, x{x—X){x — 2){x — S) .4. „ ., 

+ -^ / o o A «^ ^ ' COS*-*^ 

Wenn man die geraden Potenzen von sin t durch cos t ausdrückt; 
so findet man eine Gleichung von der Form 

cos xt = Aa: COS* t + Bg COS*"^ < + Cx COS*~* t -] 
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Es sollen die Koeffizienten Äx, B^, 0^, ».. bestimmt werden. 
Der Kürze halber wollen wir cos t mit u bezeichnen. In der 
Entwicklung 

cos xt = AxW^ + BxU*-^ + CxU*-^ + BxU'-^ + • • • 

hat die Funktion 

Ax einen Sinn fttr x^l, 

Bx }) )) >; ^ ^ 2, 

^X )y )} V ^ ^ ^9 

Bx n f, f) OS^ O, 



Aus der Formel cos t = A^u folgt, daß Ax==^ 1 ist. Um die 
übrigen Anfangswerte zu bestimmen, genügt es, in den Formeln 

cos 2t = A^u^ + JB^, 

cos 4^ = A^u^ + B^u^ + ^4, 

cos 6t = A^u^ + BßW* + CgW^ + A 



t = — oder w = zu setzen. Maji findet 

B^ = — 1, C/4 = 1 , -Dg = — 1 , . . . 

Aus der Relation 

cos (x + l)t + cos (x — 1)< = 2 cos t ' cos xt 

folgt die Identität 

Ax+iu'+'^ + (Bx+i + Ax^i)u-'^ + (a+i + Bx^i)u'-^ + . . . 

= 24rW«+^ + 2BxU"'^ + 2CxU'^''^ + • • • 
Der Vergleich der Koeffizienten fahrt zu den Differenzengleichungen 

Ax+i-2Ax^0, 
Bx+i— 2Bx= — Ax—u 
C«+i — 2 Cjc = — Bx—i, 
Bx+i--2Bx=--Cx^u 



Diese Gleichungen yereinfachen sich, wenn man neue Funktionen 
einführt. 

Wenn wir setzen 

Ax=2*-'^'ax, Ä=-2^-»-6^, a=2»-'^.c„..., 

so' finden wir ^ 

aa.4.1— aa.= U, 

&a:+l — &« = Ä;e— ;i, 
Cx+1— Cx = &ar— 1; 
da?-fl — da? = Cx—lf 
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Die Anfangswerte sind 

«1=1, &8 = 2, ^4=2, de = 2, ... 

Ans der ersten Gleichung folgte daß ax konstant ist; da aber 
«1=1 ist, so ist 

aa,= l nnd Ax=2*-^. 

Die zweite Gleichnng gibt 

Weil \^2 ist, so ist c = 0, 

hx^x und JBx = — 2*—' • rc. 

Die weiteren Bechnnngen sind 

.,= 2(.-l) = (^:4>^ + c, 

(a!-l)(ai-2) ^ 
<'»— 1.2 ^' 

n 4-8-2 /. , , 

, _ (a;-2)(a;-8)(a;-4) . .x 

"' 1.23 ^^ "" ^' 

"~ 1.2.34 1-2 "^^^ 

O 5.4.3.2 3-2, rv 

2==r72T8T4"-rT2 + ^> ^^0' 



_ (j?~3)(a;-4)(a;~5)(a?~6) (a;~5)(a;~6) 
1.2. 3. 4 12 



+ C; 



'"'"'^^L 1.2.3.4 1.2 J 

_ (a;~4) (a;-6) (a?~- 6) (a?- 7) (a;-8) (a;-6) (a?~ 7) (a? -8) 

12-3. 4. 5 123 

Q 6.6.4.3.2 4.3-2 , __ ^ 

^-1.2.3.4.6~ir2T3"*"^' ^ "" ^' 

^ ^ (a;-4)(a;-6)(a;-6)(a;-7)(a;-8) (a; - 6) (a; - 7) (a; ~ 8) 
'*"" 12. 3. 4. 5 12.3 ' 



Die gefundenen Ausdrücke für c«, (2x; ^x; /« lassen sich ver- 
einfachen. Es ist 
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«.-(i-f)(i-i)-i-^' 

rf,_(»_4).[(l-|)(l-|)-l]-2.<'-'><— " 



f 



€x — 



12-3 

{x — 6)(x — ß) r/^ ^\ /i ^\ i1 x{x — 6){x — ei)(x — l) 



[(l-f)(l -?)-'] 



1-2 L\ 3/ \ 4/ J 12.34 

(a; — 6) (rc — 7) (a; — 8) r/^ a; \ /^ a; \ ^ "] a;(aj— 6) (x—l) (x—S) {x—d) 



/■.= '"-"r7.r'" [('-!)(i-i)-']= 



1-2. 3. 4. 5 



Die gesuchten Koeffizienten haben also die Werte 

7) _-. 9x-7 a?(a?-4)(a;-6) ^ . <^ ^ a?(a; - 6) (a; - 6) (a; - 7) 
^* "^ TYTS ' ^^0.= ^' 1.2. 3. 4 ' 

jP ^_ Oa>~ii a;(a?-6)(a;-7)(a;-8)(a;-9) 
* * 1-2. 3-4. 6 

P. L. Tschebyscheff hat bemerkt^ daß die ganze Funktion 
x^ Grades 

C08 X t 

f(^) == -^r:T = w* + jp,w«-i + i>jW^-»+ . • • + p^__^u +p^ 

im Intervall (— 1, + 1) weniger von Null abweicht, als jede 
andre ganze Funktion von derselben Form: 

g)(u) = Ujc+ äiW*-^ + äaM*-^ H 1- a^-i^ + ix' 

(Oeuvres de P. L. Tchebychef, t. I, p. 295—301. Vgl. J. Ber- 
trand, Traite de calcul differentiel et de calcul integral, t. I, 
p. 512— 521.) 

Die Abweichung von Null ist der größte absolute Betrag 
einer Funktion. Wenn z. B. 

\F(u)\^Ä für a<u<b, 

so ist Ä die Abweichung von Null der Funktion F(u) im Intervall (a, 6). 
Die Abweichung von Null der Funktion 

f(u) = ^_^ cos (x arc cos u) 

im Intervall (— 1, + 1) ist gleich -^irr * 

Wir wollen voraussetzen, daß im Intervall (—1, +1) die 
Funktion g)(u) weniger von Null abweicht, als f(u). 
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Für 



ist 



' X X X 

u==l, COS — f cos — )''-> cos ^^ —f — 1, 

' X X X ' 



Wenn wir der Kürze wegen die Bezeichnung einführen 

Uk = cos — i (Ä; = 0, 1, 2, . . . , x)y 

X 

so bestehen nach nnsrer Voraussetzung die Ungleichheiten 

In jedem der Intervalle 

hat die Gleichung /*(w) — gp(M) = mindestens eine Wurzel und also 
im Intervall (—1, +1) mindestens x verschiedene Wurzeln. Dies ist 
aber unmöglich^ weil die ganze Funktion 

den Grad a? — 1 nicht übersteigen kann. 

Damit ist also der Satz von Tschebjscheff bewiesen. Diesen 
Beweis verdanke ich einer mündlichen Mitteilung von Herrn 
A. A; Markoff. • * 

(Vgl. Lacroix^ Trait^ du calcul diffi^rentiel et du calcul integral. 
Nr. 1052—1054) 

Drittes Kapitel. 
Lineare Differenzengleiehnngen mit konstanten Koeffizienten. 

54. Partikuläre LSaungen der homogenen Oloiehung. 

Um die Lösungen der Gleichung mit konstanten Koeffizienten 

ZU finden^ werden wir in die linke Seite dieser Gleichung einsetzen 

(2) y^=a^.w^, 

wo Ux eine ganze Punktion m^^ Grades ist. 
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Mit Benntznng der Formeln 

findet man ab Besnltat dieser Substitution die Identität 

+ a*+« [w(w- 1) a«-H (w-l)(n-2)i>ia«-» + (w-2)(w-3) jp^a«-* H- . .+ 2i>^_5,] ^ 

+ a*+« • w(w - 1) (w - 2) • • . 2 . 1 . jTF"^- 

Zur Abkürzung woUen wir die Bezeichnung einführen 
(3) f{d) ^^ + pi«"-^ + A^-* +••• + !>„. 

Die gefundene Identität nimmt dann die Form an 

i y,^n + Piyx^n^l + P^yx-^n^% + -' + PnVx 



X 



(4) 



= a''f(a)-u, + «•+i./*(o). A«. + o«+«/-"(o) • ^ + 



Unter v ist hier a^^u» zu verstehen. 

Die rechte Seite Ton (4) verschwindet, wenn a eine Wurzel der 
Gleidbang f{z) = und u« =" 1 ist. Eine partikulare Lösung von (1) 
ist dann a^ 

Wir wollen voraussetzen, daß a eine ifc-fache Wurzel der Gleichung 
f{z) « ist. Die rechte Seite von (4) ist dami 

A*»w 

Dieser Ausdruck ist gleich Null, sobald u» eine ganze Funktion 
vom {k — 1)*®*^ oder niedrigeren Grade ist. Wir setzen 
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Mo = 0^ AWa:=l; 



Dann ist 

Wj =1, Wa; = a?, Wa? = iz;(a? — 1), ... 

Ux = ir(ic — l)(x — 2) . . . (a; — Ä + 2); 

wir haben also h partikuläre Lösungen von (1) 

a*, a'^'X, a'^x^x— 1), . . ., a''^x{x — 1) , . . (x — Jc + 2), 

In Nr. 48 war es schon gezeigt^ daß den partikulären Lösungen 
einer linearen homogenen Gleichung ein konstanter Faktor hinzu- 
gefügt werden kann. Deswegen kann man die gefundenen partikulären 
Lösungen durch die folgenden ersetzen 

a*, ^xa'-^ "-xix — l)a*-«, . . ., ^{x — 1) . . . (x — 1 + 2)a*-*'^^ 

Um diese Lösungen zu bilden^ genügt es^ die Funktion ^ zu 
nehmen und (Je — 1) sukzessiye Ableitungen nach t zu bilden und 
in den Resultaten t durch a zu ersetzen. 

55. Allgemeine LSaung der homogenen Oleichnng. Bei 

der Auflösung der Gleichung 

werden wir zwei Falle unterscheiden: 1. die entsprechende algebraische 
Gleichung 

hat verschiedene Wurzeln und 2. die Gleichung f(0) =» hat gleiche 
Wurzeln. 

Wenn wir im ersten Falle die Wurzeln mit 

bezeichnen; dann sind nach Nr. 54 

^v ^f; • • • ; ^n 
partikuläre Lösungen von (1) und nach Nr. 48 genügt 
(5) y^- C^al + C^al+ G^al+^^^+ Cnal 

der Gleichung bei willkürlichen Werten der Konstanten 

^V ^2> ^8? • • •; ^«* 
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Um zu beweisen, daß der Ansdruck (5) die allgemeine Lösung 
ist, genügt es zu zeigen, daß Ci, Og, . . ., C« sich durch willkürlich 
gegebene Werte y^, y^, y^^ . . ., y^_^ ausdrücken lassen. 

Um das System 



(6) 



.y«-i =- Ci«r '+ c'aar '+ Cs«r '+ • • • + (^nai-^ 



aufzulösen, werden wir die Gleichungen des Systems mit unbestimmten 
Paktoren g„__i, Q-n^v ^n—v • • • ? ^0 multiplizieren und dann addieren. 
Das Resultat nimmt die Form an 

wenn man die Bezeichnung 

(8) (p(0) = go^~'+ 3i ^""* + • + q„^^0 + q^^^ 

einführt. 

Um Cj zu bestimmen, werden wir g^, ffi, . . . q^—i ®^ wählen, daß 

9^K) = 1^ 9K) == 0; 9(%) = 0, . . . , y (an) = 
wird. Es ist in diesem Falle 

''^^^~" («i-a2)(«i-as)... («!-«„) 

und die Relation (7) gibt 

In derselben Weise lassen sich C^, Cg, ..., (7„ durch y^, y^, 
y^y • • •> ^n-i ausdrücken. 

Damit ist also bewiesen, daß der Ausdruck (5) wirklich all- 
gemeine Lösung von (1) ist. 

Jetzt wollen wir voraussetzen, daß 

gleiche Wurzehi besitzt. 

Li allen einzehien J^ällen ist die Beweisführung dieselbe. Des- 
wegen genügt es, nur einen Fall zu betrachten. Wir werden voraus- 
setzen, daß 

f(0) = (^ - ai)\z- a^ijs - ag) . . . (^ - a^ 



(9) 



(10) 
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ist und daß die Zahlen a^, a^y c^, .,., an voneinander ver- 
schieden sind. 

Da wir n partikuläre Lösungen von (1) kennen (Nr. 54), so 
genügt der Ausdruck 

unsrer Differenzengleichung. Es ist noch zu zeigen, daß das System 

y^^C^a^ +C2-1 +C3O -^-C^a^ + +0„a„, 

yo-C^'l +C^'0 +C3.O -\-C^l + +C„.l 

in bezug auf C^, Cg, . . ., 0„ lösbar ist. 

Wenn man die Gleichungen (10) mit g^, ^i, ..., Q^—v 9n—i 
multipliziert und zusammenaddiert, so findet man 

1= Ci<p(a,) + 02<p'(aO + C3<)p''(a0 + C,<)p(a4) + . • . + 
Hier ist 

9 W = 2o^"~^ + «1^-^ + • + g„_2^ + g„__r 

Um Ci zu bestimmen, genügt es, die Koeffizienten g^, q^, . . ., 
g„_2' ««-1 ^^ 2u wäMen, daß 
VK) ^ 1> ^^'K) = 0, y"(ai) = 0, ^(aj = 0, . . . , q>{an) = 0. 
Um die Gleichungen 

zu befriedigen, genügt es, zu setzen 

(12) (p(0) = (^ - aj(^ - «5) • • • (^ - a«)^(.^); 

wo ^(;8f) eine willkürliche Funktion ist. 
Aus (12) und aus den Bedingungen 

(|p(a,) = l, ()p'(aO = 0, ()p"(ai) = 

folgen die Werte für ^(a^), ^'(a^), ^"(«i). Die Funktion vom 
niedrigsten Grade, die allen diesen Forderungen genügt, ist 

y(^) = (^ - aj(^ - Og) . . . (^ - a«)[^(ai) + (^ - ai)^'(ai) 

Seliwanoff, Düferenzenreohnung. 6 
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In derselben Weise lassen sich die Eonstanten O^ und C, durch 
Vof ViJ "'f y„-i ausdrücken. 

Um 0^ zu bestimmen, genügt es, zu setzen 



vW 



(g~a^)»(g-aa)(g~ae)...(g-a„) 
(«4 - öi)'(«4 - «ö) («4 - «e) • • • («4- «n) 



Ebenso lassen sich die übrigen Eonstanten C^, C^, ..,, Cn 
bestimmen. 

Damit ist also bewiesen, daß (9) wirklich eine allgemeine 
Lösung der vorgelegten Differenzengleichung ist. 

Durch die Einführung andrer willkürlicher Eonstanten kann 
diese Lösung eine einfachere Form annehmen. Da der Ausdruck 

Ol + C^a-^x + Csa-^x(x - 1) 

eine willkürliche Funktion 2*^^^ Grrades ist, so kann man ihn durch 
J^ + A^x + ÄqX^ ersetzen. 

Wir bekommen die folgende allgemeine Lösung unsrer Gleichung 



(13) 



y^=(A^+A^x + J«ic»)aj + C^aJ + • • • + C^a'^. 



In allen Fällen gibt es also solche n Lösungen u^^\ u^*\ . . ., w^"' 
der gegebenen homogenen Gleichung, daß man das System 



y„_, = G.u^U + <7,«« 1 + ■ ■ • + (?„«^ii 



in bezug auf C^, C^, ..., 
Determinante 



(14) 



<) 
«« 



Cn auflösen kann. Deswegen ist die 



«<*) ... «W 
mW ... M<") 



mW mW .... m^) , 



von Null verschieden. 

Als Beispiele wollen wir einige Gleichungen mit ihren allgemeinen 
Lösungen angeben. 

y= Gl -2*+ (7g -3- 
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56. G-leichnngen mit reellen Xoefflzienten. Wenn alle 
Koeffizienten der Gleichung 

reell sind^ so kann man verlangen; daß in der allgemeinen Lösung 
keine imaginären Ghrößen auftreten. 

Wir werden Yoraussetzen^ daß die entsprecliende algebraische 
Gleichung 

eine dreifache komplexe Wurzel hat^ die wir mit a^ bezeichnen. 
Nach einem Satze der Algebra hat diese Gleichung eine dreifache 
konjugierte Wurzel, die wir mit a^ bezeichnen. Die übrigen Wurzeln 
d^y a^f , > ., an sollen reell und voneinander verschieden sein. 
Die allgemeine Lösung unsrer Gleichung ist 

Wenn 

a^ = p (cos « + i sin a) 
gesetzt wird, so ist 

«4=9 (cos a — i sin a), 

aj = ^* (cos xa + i sin xa), 

aj = p* (cos xa — ■ i sin xa). 
Deswegen ist 

(Ci + C^x + C.x'^af + {C^+ C,x + G^x')al 
- [(Ci + ^2^ + Cs^ +(C^+C^x+ C.x'^^Q^ cos ax 
+ [(Ci + C^x + C^x^i - (0^+ C^x + Cf^s^)q ^ sin ax. 
Wenn man neue Eonstanten einfährt 

so erscheint die allgemeine Lösung unsrer Gleichung in der Form 
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Vx"^^' M + A^ + -4^*) <5os xa + {B + B^x + B^x^) sin xa\ 

+ C,a? + .- + CX, 

die keine imaginären Größen eniMlt. 
Beispiele: 

y. = 2» . [^ coB (a; ¥) + B sin {x ^)]- 

y,= Gl cos (a; jj + C, sin^a;^) + CjCob (a;^) + Q 8in(a;^V 

57. ToUatfindlge Oleiehnng. Um die Losnng der toU- 
standigen Gleichung 

(15) y,+„+i'iy«+„_i+i>,y,+„_,+ • • +i'„y«= Q, 

zu suchen, werden wir Yoraassetzen, daß die Losung der entsprechenden 
homogenen Gleichung 

schon bekannt ist. Es sei die allgemeine Lösung dieser Gleichung 

Statt der willkürlichen Eonstanten werden wir unbekannte 
Funktionen von x einführen^ und der Ausdruck 

(16) y = (7(1)^(1) + (7(^>mW + . . . + (7(-)w(»> 

'^ ' *^ X X X X X XX 

soll der vollständigen Gleichung genügen. Eine Bedingung haben 
wir zu erfüllen und n Funktionen zur Verfügung. Deswegen können 
wir n — 1 neue Bedingungen einführen. Lagrange wählt dieselben 
so, daß y^i^f VxA-i' ' ' '' ^ar+n— 1 ^^i^selbe Form haben sollen, als ob 
die Cx Konstante wären. 

Da Cx+i = Cx + ^Cx ist, so finden wir das System von Gleichungen 

t/ = n(i) . w(i) 4- c(*) . WW -^ f- C^*») • uH 
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Diese Formeln gelten^ wenn die Funktionen Gx den Bedingungen 
genügen 



(17) 



^%, --^^i^'+w^V» •Acrw+-.- + t*^;, .A0(«)=0, 



«<!) ,-ACW + «<« .ACW+--+»2'1- ,-AC'W=0. 
Es ist 

+ «»). „ • ACW + u^n . ACW + • • • + «ä.« • ACW, 

• x-j;-n X ' ay-f-w x " ' aj-f-n a; ' 

•v 

und es kann nicht mehr y , dieselbe Form haben, als ob die Gx 
konstant wären. 

Wenn wir die Werte y^, y^^^, y^^^^ . . ., y^^„ in die gegebene 

vollständige Gleichung einsetzen, so fallen die Glieder mit G^^y 
G^\ . . . , G^^ weg und es bleibt die Gleichung 

(18) t*u^„ . Aqy + «0.)^^ . Aq«) + . . . + «w ^ . ACW = «,. 

Aus den n Gleichungen (17) und (18) mit den n Unbekannten AC^^), 
^Gjp, . . . ; AO^") lassen sich dieselben bestimmen, wenn die Determinante 



F(x) = 









mW. 

mW. 



«1% €4« • • • «?;» 



von Null verschieden ist. 

Wir haben schon gesehen (Nr. 55), daß die partikulären Losungen 
^« ^ ^T> ' ' -> ^L**^ ^®^ homogenen Gleichung so gewählt werden können, 
daß die Determinante 



Fi-l) 



<) «<«) ... «w 

mW mW ... mW 






nicht gleich Null ist. 
Aas der Relation 

folgt, daß 
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Fiz)^ 









«Iris 



"ar+n— 1 



x+n — 1 



z-^n — 1 



— P^ ^« ^ — jp« w?^ • • • — pM^^ 

weil man za den Elementen der letzten Zeile die Elemente der vorigen 
Zeilen, resp. multipliziert mit p^, p^, . . ., p^_^y addieren kann. 
Es besteht also die Relation 

F^x)^±p^.F{x^l). 

Wäre F{x) = 0, so müßte auch F(x — 1) = sein, weil p^ nicht 
NuU ist. 

Aus denselben Ghründen wäre auch 

was aber nicht der Fall ist. 

Es ist also bewiesen, daß F(x) nicht Null ist. Die Gleichungen 
(17) und (18) sind lösbar. Man findet 

X X ' X X ' ' X X * 

Die Summation gibt Resultate von der Form 

(7(1)= «;(!)+ (7^, C(») = «4»)+C„..., q-)=«4«)+o„. 

Demnach ist die Lösung der vollständigen Gleichung 

+ w(i) «;(!) + m(2) tc;(8) ^ 1- mW t(;(~). 

X a; X X 'XX 

Es ist die allgemeine Lösung, weil die Determinante (14) von 
Null verschieden ist. 

Diese Methode von Lagrange heißt Variation der will- 
kürlichen Konstanten. (Lagrange, Oeuvres, t. 4, p. 156. Recherches 
sur les series recurrentes.) 

Als Beispiel nehmen wir die Gleichung 



Wir setzen 

Es ist 



wenn 



yxH-2-7y^+i+6yx=^- 
ACW+6'+^-A(7W=0 



ist. 
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Es wird y^ der Gleichung genügen^ wenn außerdem 

ist. Daraus folgt 

Die Summation gibt 

Die allgemeine Lösung unsrer Gleichung ist also 

58. Man kann noch in einer andern Weise von der Lösung 
der homogenen Gleichung zur Yollständigen (15) übergehen. Es 
genügt, eine partiknlare Lösung der yollständigen Gleichung zu 
kennen (Nr. 48). Wenn das letzte Glied 

ist; wo Vx eine ganze Funktion ist, dann ist es leicht^ eine solche 
Lösung zu finden. 

Wir setzen — 

wo Ux eine ganze Funktion ist; und führen diese Größe statt y^ in 
die vollständige Gleichung ein. Vermöge der Identität (4) findet 
man die Gleichung 

. . . + a'^-n . ^(«)(^) . j-^^ = a' . Vx, 

die zur Bestimmung von Ux dient. 

Der Faktor a* kann gehoben werden und es wird diese 
Gleichung die Form annehmen 

A*w A**w 

(19) /•(a)M,+a/-'(a)A«,+aY''(o)-Y72 +-+«"/'^"W TTä::^»«.- 
Wenn a keine Wurzel der Gleichung 

ist, dann ist Ux von demselben Grade wie v«. 
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Es sei 

___ a; (gg — 1) (jg — 2) . . . (a; — m -f 1) . a;(a; ~ 1) . . . (a? — w -f 2) 
^'~ ^ 1.2.3...m ^ ^ 1.2.3...(w-l) 

, ir(a;— l)...(aj — w + 3) , , x{x — l) , , 

+ "« l-2...(m-2) +--- + «».-«-IT2- + «"-^^ + '^> 

„ a;(ae — !)(« — 2).. .{x — m-\-\) , » «(as — 1) . . . (« — t» + 2) 

^*~ P» I.2.8...»» •" ''i 1.2-3...(»-l) 

+ P« l-2...(w-2) + • • • + Pm-2 1.3 + Pm-ia; + /»m- 

Die Zahlen /3 sind gegeben nnd die a werden gesnclit. 
Ans (19) findet man durch den Vergleich der Koeffizienten 

«o/'(«) =» ßt>, 
fhfifl) + ao«/"'(ö) = A» 

«,/-(a) + ««,ar («) + a,o«^ + «««'S^ = ft 



Da f{d) nicht Null ist^ so bestimmen sich sukzessive die Zahlen 

^o> ^> ^; ^> • • • 

Wenn a eine ä; fache Wurzel von f{z) = ist, so nimmt die 

Grleichung (19) die Form an 

Die linke Seite dieser Gleichung ist eine ganze Funktion 
^ten Grades, wenn Ux vom (Je + w)*^ Grade ist. 
Man kann setzen 

05 (aj — 1) . . . (a; — & — m + 1) , a; (a; — 1) . . . (a; — Ä — w + 2) , 

^'~*^ 1.2.3.. .(Ä + w) *■ **i 1.2.3...(A; + m-l) *" 



» • • 



a;(a; — 1). . .(a: — Ä;) , aj(a; — 1) . . . Ca? — Ä + 1) 
"^ ^-^ 1.2.3...(* + 1) "*" ^ 1 . 2 . 3 . . . Ä; ' 

und zur Bestimmung der Koeffizienten findet man die Gleichungen 

"«" l-2...*~Po' 

'^" 1.2... * + **!* 1.2...(ik+l) + *^* 1.2...(Ä + 2) «• 
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Diese Gleicliimgen sind lösbar ^ weil a und f^^^(a) nicht gleich 
Null sind. 

Z. B. wollen wir eine partikulare Lösung der Gleichung 

suchen. Hier ist a == 1 und dieser Wert ist eine einfache Wurzel 
der Gleichung 

Deswegen gibt es eine partikuläre Lösung von der Form 

Li die gegebene Gleichung setzen wir ein 

und wir findeiju % 

,3 1 

*^6Ö' ^ = -10- 

Nach dem bekannten Satze (Nr. 48) ist die allgemeine Lösung 
dieser Gleichung 

Die Gleichung 
hat eine partikulare Lösung von der Form 

Xx = äq • 

Aus der Gleichung 

folgt; daß ckq =» ö^ ist. Deswegen ist die allgemeine Lösung unsrer 
Gleichung 

Wir wollen jetzt einige Anwendungen der Differenzengleichungen 
höherer Ordnung zeigen. 

59. Bekurrente Beihen. Dies sind solche Reihen 

« « ... 

Vq) Vi) y%y • • •; y«; ^x+i' * * *' 

für die je n aufeinander folgende Glieder einer Relation 
mit konstanten Koeffizienten genügen. 

Seliwanoff, Differenzenreoluitixig. 6** 
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Die Anflösnng dieser DifiPerenzengleichimg dient dazu^ mu das 
aUgemeine Glied y^ als Funktion von x darzustellen. 
Als Beispiel wollen wir die Reihe betrachten 

0, 1, 1, 2, 3, 5, 8,..., 

deren jedes Glied die Summe der beiden vorhergehenden Glieder ist. Je 
drei nacheinander folgende Glieder sind durch die Relation verbunden 

Die Anfangswerte sind 

Die Wurzeln der quadratischen Gleichung 

sind ^— ^— 

l±f. und 1^. 

Deswegen ist die allgemeine Lösung der Differenzengleichnng 

Wenn C^ und C^ durch die Anfangswerte Pq und y^ bestimmt 
werden; findet man^ daß das allgemeine Glied unsrer Reihe 

ist. 

Die Potenzreihe 

F(t)=^yo+ytt + y,t^+y,t' + -' + yJ' + ^'', 

deren Koeffizienten eine rekurrente Reihe bilden^ ist gleich einer 
rationalen Funktion von & Die Relation zwischen den n Koeffizienten 
der Reihe sei 

Nach Multiplikation von F(t) mit 
findet man 

Vo + {yi+Piyo)^+ (j/a +Piyi+P2yo)^^ + (ys +Piyi H-jp^yi +PByo)^^ + 

Die Koeffizienten der folgenden Potenzen sind die Werte der 
Funktion 

für x = 0, 1, 2; . . . Deswegen sind alle diese Koeffizienten gleich Null. 



Drittes Kapitel. Lineare Differenzengleichtmgen mit konstanten Koeffizienten. 91 
Umgekehrt, die rationale Funktion 

läßt sich in die Beihe 

entwickeln, für die je n aufeinander folgende Koeffizienten durch die 
Relation 

verbunden sind. Die n ersten Koeffizienten bestimmen sich aus den 
Gleichungen 



Also können die Differenzengleichungen dazu dienen, um die 
Koeffizienten der Entwicklungen rationaler Fimktionen in Potenz- 
reihen zu bestimmen. 



60. TschebyscheflQeicher Beweis eines Saties von Lfun6. 

In seinen Vorlesungen zeigte Tschebyscheff die folgende An- 
wendung der Differenzengleichungen. 

Man soll bestimmen, wie groß die Anzahl der sukzessiven 
Divisionen sein kann, um den größten gemeinschaftlichen Teiler 
zweier Zahlen a und h zu finden. 

Es sei 

Die Anzahl der Divisionen ist hier n. 

Weil alle Quotienten größer oder gleich Eins sind, so ist 

. . . Tn—Z ^ ^n— 2 + **«—!> ^n-2 ^ ^n — 1; ^n— 1 ^ 1- 

Daraus folgen die Beziehungen 

r„=0 r„_i^l, r„_2^1, rn_s ^ 2, 

^n— 4 ^ 3, rn— 5 ^ 5; ^n— 6 ^8, ... 
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Wenn wir die Zahlen 

2/0 = 0, ^1=1, % = !, ya, = 2, ..., ya.= y^^i + y^_a, ... 

einführen, so ist 

Wenn für ein^. ganze positive Zahl m die Ungleichheit 

besteht, dann ist notwendig y^ > y^ und m > w. Die gesuchte 
Anzahl der Divisionen ist also kleiner als eine solche Zahl m, für die 

ist. Den Ausdruck für v kennen wir schon aus Nr. 59. Es ist i 

•C I 

deswegen die Ungleichheit j 

aufzulösen. 
Da 

|/5==2,2...; Ii5^.==l, 6...; ill»^ = - 0, 6 . . ., 
so können wir y^ annähernd in der Weise ausdrücken 

wo 9 zwischen und 1 ist. 
Wenn 

ist, dann ist 

y„^6 + (l±0) und folglich y^>6. 
Es genügt 

^ log [2, 2. (5 + 1)] 

zu nehmen. Wenn die Zahl b aus k Ziffern besteht, dann ist der 
gemeine Logarithmus von 2, 2-(& + 1) ungefähr gleich Ä, log 1, 6 

ist beinahe — und 

^ m^5k. 

Daraus folgt der Satz von G.Lame! Die Anzahl der Divi- 
sionen, die auszuführen ist, um den größten gemeinschaft- 
lichen Teiler der Zahlen a und & (a> 6) zu finden,- ist kleiner 
als die fünffache Anzahl der Ziffern der Zahl b. 

Lame hat diesen Satz auf elementarem Wege bewiesen. (Paris. 
Comptes rendus, 19. [1844], p. 867.) 
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Professor in Strafibnrg Professor in Oiefien. 

In 3 Bänden. Band I (XTV n. 446 S.) gr. 8. 1903. In Leinw. geb. n. •# 8.— 

Das Werk, dessen erster Band soeben erschienen ist, richtet sich in erster 
Linie an die Lehrer, die darin Anregung finden sollen, ihren Unterrichtsstoff 
auszuwählen und namentlich in den höheren Klassen zu vertiefen, sodann aber 
auch an Studierende, die eine Anlehnung an die Elemente und Auffrischung und 
Ergänzung früher erworbener Kenntnisse suchen. 

Durch das Zusammenwirken mehrerer Grelehrter hoffen die Herausgeber, die 
möglichste Vollständigkeit zu erreichen. Der erste Band umfaßt den algebraisch- 
• analytischen Teil. Der zweite Band, der Ostern 1904 erscheint, wird die Geometrie 
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entnommen ist. Die Vorarbeiten sind so weit gediehen, daß die Vollendung des 
ganzen Werkes im nächsten Jahre erwartet werden darf. 



